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Abstract. We investigate the properties of the extinction curve in the rest-frame UV for a sample of 34 UV-luminous galaxies 
at 2 < z < 2.5, selected from the FORS Deep Field (FDF) spectroscopic survey. A new parametric description of the rest- 
frame UV spectral energy distribution is adopted; its sensitivity to properties of the stellar populations or of dust attenuation 
is established with the use of models. The latter are computed by combining composite stellar population models and calcu- 
lations of radiative transfer of the stellar and scattered radiation through the dusty interstellar medium (ISM) for a dust/stars 
configuration describing dust attenuation in local starbursts. In the favoured configuration the stars are enveloped by a shell 
with a two-phase, clumpy, dusty ISM. The distribution of the z ~ 2 UV-luminous FDF galaxies in several diagnostic diagrams 
shows that their extinction curves range between those typical of the Small and Large Magellanic Clouds (SMC and LMC, 
respectively). For the majority of strongly reddened objects having a UV continuum slope /? > -0.4 a significant 2175 A ab- 
sorption feature (or "UV bump") is inferred, indicating an LMC-like extinction curve. On the other hand, the UV continua of 
the least reddened objects are mostly consistent with SMC-like extinction curves, lacking a significant UV bump, as for the 
sample of local starbursts investigated by Calzetti and collaborators. Furthermore, the most opaque (fi ~ 0) and, thus (for our 
models), dustiest UV-luminous FDF galaxies tend to be among the most metal rich, most massive, and largest systems at z ~ 2, 
indicating (Z) ~ 0.5 - 1 Zq, (Mstars) ~ 6 x 10'" Mq, and (Rcn) ~ 4kpc, respectively. The presence of the UV bump does not 
seem to depend on the total metallicity, as given by the equivalent width (EW) of the C IV doublet. Conversely, it seems to be 
associated with a large average EW of the six most prominent interstellar low-ionisation absorption lines falling in the FORS 
spectra. The average EW of these saturated lines offers a proxy for the ISM topology. We interpret these results as the evidence 
for a difference in the properties of the dusty ISM among the most evolved UV-luminous, massive galaxies at z ~ 2. 

Key words, galaxies: high-redshift ~ galaxies: starburst - galaxies: ISM - dust, extinction - ultraviolet: galaxies 



J> 1 . Introduction The direct determination of the spectral dependence of ex- 
tinction from observational data is limited to our own Galaxy, 

The electromagnetic radiation propagating through a medium the Magellanic Clouds, and M 3 1 (for reviews see Fitzpatrick 

containing dust particles experiences extinction. In general, the 120041 and Clavton l2004> . 

intensity of the transmitted beam is reduced by two physical c t /-^ i f ^ a 

■' J f J Extinction curves or our own Galaxy vary for difrerent 

processes: absorption and scattering. Dust grains absorb and • u» r j * j j ^ » • » /f j 

^ ^ & 6 sight lines and towards difrerent environments (from dense 

scatter more the shorter the wavelength. As a consequence, , , i j ^ ,.u j fr * n a- 

° ^ molecular clouds to the diftuse interstellar medium, ISM; see 

both the total intensity and the spectrum of the pristine radi- . • i Hnn/i t u j » • j 

J f Fitzpatrick '20041). However, the extinction curves determined 

ation are modified. The so-called extinction curve quantifies ^ xc,.. ^ ^, ,•^mT^ 

^ for the diffuse ISM of the Milky Way (MW) seem to repre- 

the amount of extinction as a function of wavelength along the , . f i . u • * . i . 

° ° sent a one-parameter family, the parameter being the total-to- 
line of sight. It depends on the chemical composition, structure, , ^. ^. n A i-n rr^ A ^ loon 
° ^ 1- ' ' selective extinction ratio Kv = Av/£,g_y (Cardelli et al. 1989, 

^fr^tl distribution of the in™^ dust grains (Mathis et |^ Fitzpatrick p99 , Clayton et al. m Valencic et al. 

al. n977.. Desert et al. 1990 Dweck et al. 1997. Draine & Li . . • i onn/i > d ^7"^ i- v ■ 

, — ' ' 2004 , see however Fitzpatrick 2004 . Broadly speaking, Ry is a 

EUnn Li & Draine 2001 , Zubko et al. 2004. Hence, studies of . ■ ■ ^ „ ■ Lv u 

' ' ' — ' proxy for the average grain size: low-Ky sight lines have more 

extmction provide information on the properties of dust. „ • »u u- u d • u» r t-u n »u 

^ ^ ^ small grains than mgn-Rv sight lines. The smaller the grains, 

the steeper the extinction curve in the ultraviolet (UV). 



Send offprint requests to: S. Noll The average extinction curve for the diffuse ISM of the 

* Based on observations obtained with FORS at the VLT, Paranal, MW has /^i/ = 3.1. It is characterised by the presence of a 

Chile on the course of the observing proposals 63.0-0005,64.0-0149, bump at a constant central w avelen gth of 2175 A. This ab- 

64.O-0158, 65.O-0049, 66.A-0547, 68.A-0013, 68.A-0014, 69.A- sorption feature (Witt & Lillie'1973) is most likely produced 

0104 by graphitic carbon, though its origin is still under discussion 
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(DraineEOnil Whittet lMni and references therein). The width 
of the UV bump seems to depend on the environment: narrow 
bumps tend to be associated with diffuse regions, while broad 
bumps favour dense regions (Valencic et al. 2004 1. The occur- 
rence of weak bumps may reflect processing that modifies or 
destroys the bump carrier in dark clouds (Whittet et al. 2004). 

The few sight lines studied in M 3 1 seem to show a CCMl 
far-UV extinction and a weak 2175 A bump (Bianchi et al. 
I1996> . There are remarkable differences in dust grain properties 
between the Milky Way and the Large and Small Magellanic 
Clouds (LMC and SMC, respectively). The sight lines toward 
the LMC 2 supershell and the average LMC extinction curve 
exhibit a weaker bump and stronger far-UV extinction than the 
average MW Ry - 3.1 extinction curve. The sight lines toward 
the SMC bar exhibit extremely weak or absent UV-bump ab- 
sorption and a very strong far-UV extinction, which is indica- 
tive of a very large fraction of small dust particles (Rodrigues et 
al.[l997 1. These diff'erences may reflect differences in metallic- 
ity' and gas-to-dust ratio (Bohlin et al. il978i Luck & Lambert 
n^92l Gordon et al. 2003 1. 

Different environments, such as star-formation regions 
where large amounts of UV radiation and shocks are present, 
may play a fundamental role in processing dust. This is not 
easy to establish however. As remarked by Clayton (2004), the 
LMC 2 supershell is located in the 30 Dor region, which is a 
much larger star-forming region than any in the SMC, but the 
LMC 2 supershell extinction curve is in much better agreement 
with a ICCMI MW extinction curve than the SMC extinction 
curve. Nevertheless, it is intriguing that nearby starburst galax- 
ies spanning a large range in metallicity (from sub- to super- 
solar) possibly contain dust with an extinction curve lacking 
a 2175 A bump, like the SMC curve, and a steep far-UV rise, 
intermediate between the MW and SMC curves (Calzetti et al. 
[T994II2ODO1, Gordon et al. [T997l Leitherer et al. 2002). For these 
objects, the properties of dust (i.e. the extinction curve) can be 
determined only in an indirect and rough way, after resorting to 
models of radiative transfer of the stellar and scattered radiation 
through a dusty ISM (Witt &. Gordon 2000, WGOO Misselt et 
al.gOOl Pierini et al. l2004l Inoue|2005). In fact, on a galaxy 
scale, the propagation of electromagnetic radiation through a 
dusty medium is described by the attenuation function, i.e. the 
combination of the extinction curve with the geometry of the 
system, in which a substantial fraction of the scattered light is 
returned to the observer's line of sight. 

SMC-like dust seems to characterise also starburst galaxies 
in the distant universe. This holds for individual objects with 
a well-sampled SED (e.g. Gordon et al. 1999) as well as for 
large photometric samples, like that of ~ 1000 Lyman-break 
galaxies (LBGs) at z ~ 3 observed by Steidel et al. (2003) 
(Vijh et al. 2003). The absence of a 2175 A bump seems to 
characterise the rest-frame UV spectrum of other objects ly- 
ing at intermediate/high redshift, like quasars/AGN (Pitman et 
al.pOOOi Maiolino et al. 2001 Hopkins et al. 2004) and GRB- 
host galaxies (e.g. GRB 020813 at z = 1.255, SavagUo & Fall 

' The LMC and SMC are deficient in heavy elements by factors 
of about 2.5 and 7, respectively, compared with the solar standard 
(Westerlund i 19971 . 



I2004> . On the other hand, a UV bump seems to be detected 
in the composite absorption spectrum of 96 intervening Mg II 
systems at 0.2 < z < 2.2 (Malhotra 1997). However, Khare 
et al. ( I2005> come to an opposite conclusion for the compos- 
ite absorption spectrum of 815 intervening QSO absorbers at 

1 $ z ^ 1.9. Further evidence for MW-type dust being also 
present at high redshift comes from a composite spectrum of 
nine radio galaxies at z ~ 2.5 (Vernet et al. I2001> . and from 
the individual absorption spectra of three intervening Mg II ab- 
sorbers at 1.4 ^ z ^ L5 (Wang et al. .2004j . 

It is not clear when star-forming galaxies at high redshift 
exhibit an extinction curve with a 2175 A bump. Therefore, 
we have selected a sample of 34 UV-luminous galaxies at 

2 < z < 2.5 from the FORS Deep Field spectroscopic sur- 
vey (Noll et al. 2004, see Sect. |2j. We will characterise their 
composite extinction curve at rest-frame UV wavelengths in 
terms of far-UV slope and presence/absence of the UV bump. 
The UV-luminous galaxies at z ~ 2 exhibit large values of UV 
continuum reddening (e.g. van Dokkum et al. 120031 Daddi et 
al. l2004l Noll et al.f2004), if the slopes of their rest-frame UV 
continuum spectra are interpreted as for local starburst galaxies 
(Meurer et al. .l999. Calzetti et al. 2000 Leitherer et al. .20{nt . 
In this case, one can obtain attenuation factors at 1500 A up to 
~ 10^ (Noll et al. l2004t . that are more than one order of mag- 
nitude larger than those determined for nearby starbursts, and 
correspond to the largest values estimated for LBGs at z ~ 3 
by Vijh et al. (2003 ). At variance with those LBGs at z ~ 3, the 
UV-luminous galaxies at z ~ 2 are characterised by half-light 
radii of several kpc, total stellar masses between a few 10'° and 
more than 10" M©, metallicity between the LMC one and so- 
lar, and the absence of significant Lya emission (Mehlert et al. 
2002 Daddi et al.l2004l Noll et al. l2004l S hapley et al .l2004l 
Steidel et al. l2004l van Dokkum et al. l2004l Drory et al. l2005l . 
Hence, the properties of the mixture of dust grains present in 
UV-luminous galaxies at z ~ 2 may differ from those of LBGs 
at z ~ 3 and those of local starburst galaxies. 

In the following we illustrate this by making use of a new 
parametric description of the rest-frame UV SED that allows 
us to define dust-sensitive diagnostic diagrams. We interpret 
the distribution of the FDF sample galaxies in each of these 
diagrams through a comparison with that obtained either for a 
local reference sample of starburst galaxies or for models of 
starburst galaxies combining stellar population synthesis and 
radiative transfer calculations (Sect.|3}- 

Throughout this paper - 70kms"' Mpc"', Q^. = 0.7, 
and Qm = 0.3 are adopted. 

2. The spectroscopic sample: selection and data 

2.1. The FORS sample of UV-luminous galaxies at 

We have selected 34 galaxies (see Table [TJ with /vega 5 24.5 
in the spectroscopic -redshift range 2 < z < 2.5 from the cat- 
alogue of the FORS Deep Field (FDF) spectroscopic survey 
(Noll et al. l2004t . The redshift interval chosen guarantees the 
covering of the rest-frame UV wavelength range from Lya to 
around 2500 A, including the 2175 A feature, with the avail- 
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Table 1. The spectroscopic sample of FDF galaxies at 2 < z < 
2.5 (Noll et al. l2004> . IDs, world coordinates, and total / mag- 
nitudes were taken from Heidt et al. (2003 1. Redshifts were de- 
rived by Noll et al. 82004 1. Luminosities at 1500 A were taken 
from the same source, though corrected for the modified cos- 
mology used. 



No. 




RA 


DEC 




z 


log Li500 


(FDF) 


(J2000) 


(J2000) 


[mag] 




[W/A] 


1208 


01 


05 51.9 


-25 48 04 


23.68 


2.178 


33.88 


1691 


01 


05 53.7 


-25 45 29 


23.89 


2.344 


33.94 


1744 


01 


05 53.9 


-25 46 06 


24.10 


2.374 


33.91 


1991 


01 


05 54.8 


-25 46 16 


24.33 


2.334 


33.85 


2274 


01 


05 55.9 


-25 44 34 


23.34 


2.253 


34.18 


2418 


01 


05 56.4 


-25 45 12 


23.16 


2.332 


34.25 


2495 


01 


05 56.7 


-25 43 44 


23.31 


2.453 


34.25 


2636 


01 


05 57.3 


-25 44 02 


23.43 


2.253 


34.01 


3005 


01 


05 58.6 


-25 48 14 


23.51 


2.253 


34.03 


3163 


01 


05 59.2 


-25 45 38 


23.35 


2.441 


34.21 


3300 


01 


05 59.6 


-25 46 30 


23.91 


2.375 


34.07 


3374 


01 


05 59.9 


-25 45 1 1 


23.34 


2.386 


34.27 


3688 


01 


06 00.9 


-25 47 05 


24.07 


2.375 


34.00 


3810 


01 


06 01.3 


-25 45 28 


22.67 


2.372 


34.44 


3874 


01 


06 01.5 


-25 45 46 


23.30 


2.483 


34.23 


3875 


01 


06 01.5 


-25 47 34 


24.53 


2.243 


33.84 


3958 


01 


06 01.8 


-25 44 29 


23.87 


2.130 


33.97 


4795 


01 


06 04.8 


-25 47 14 


23.31 


2.159 


34.13 


4871 


01 


UO 05 . 1 


-25 46 04 


23.39 


2.472 


34.26 


4996 


01 


06 05.5 


-25 46 28 


23.25 


2.028 


34.04 


5058 


01 


06 05.7 


-25 46 26 


23.34 


2.027 


34.15 


5135 


01 


06 06.0 


-25 44 44 


23.62 


2.346 


34.02 


5165 


01 


06 06.1 


-25 44 43 


23.26 


2.346 


34.34 


5190 


01 


06 06.1 


-25 44 43 


24.39 


2.347 


33.83 


5227 


01 


06 06.3 


-25 43 52 


23.85 


2.399 


34.08 


6024 


01 


06 09.2 


-25 48 14 


22.00 


2.372 


34.68 


6372 


01 


06 10.5 


-25 48 29 


23.38 


2.349 


34.21 


6407 


01 


06 10.6 


-25 45 32 


23.59 


2.162 


33.95 


6934 


01 


06 124 


-25 44 57 


22.90 


2.445 


34.43 


6947 


01 


06 124 


-25 48 15 


23.83 


2.357 


33.97 


7029 


01 


06 12.7 


-25 45 58 


23.63 


2.374 


34.12 


7078 


01 


06 12.8 


-25 46 01 


23.97 


2.378 


34.00 


7307 


01 


06 13.6 


-25 47 25 


24.07 


2.438 


34.01 


7342 


01 


06 13.7 


-25 46 13 


23.80 


2.375 


34.06 



able optical spectroscopy. The limiting /-band magnitude pre- 
vents the contamination from faint Lyor emitters that were dis- 
covered serendipitously. With an average luminosity at 1500 A 
(rest frame) <Li5oo) = 1.4 x lO-'^^W/A (not corrected for dust 
attenuation), the present sample of 34 galaxies represents the 
bright end of the distribution in L1500 of FDF galaxies at z ~ 2 
(see Fig. Hereafter it is referred to as the FDF sample of 
"z ~ 2 UV-luminous galaxies". 

The VLT-FORS spectra of our sample galaxies were ob- 
tained with the grism 1 50 1 and a 1" slit width (Noll et al.E004 1, 
resulting in a measured spectral resolution element between 18 
and 23 A that coiTesponds to rest-frame 5 to 8 A. The spectra 
cover the entire optical wavelength range providing reasonable 
data between 3500 and 9300 A (observed frame), which allows 
us to probe the UV SED between 1216 A (i.e. the Lya transi- 
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Redshift 

Fig. 1. The luminosity at 1500 A (rest frame) L1500 as a func- 
tion of redshift for the 118 galaxies at z > 1.35 in the FDF 
spectroscopic sample (from Noll et al. "20041. The 34 galaxies 
under study (see Table0 are reproduced with circles. Among 
the remaining high-redshift objects, star-forming galaxies are 
marked by crosses and QSOs by asterisks. L1500 refers to the 
measured emission within the rest-frame wavelength range 
from 1480 to 1520 A. 



tion of hydrogen) and ~ 2500 A, bracketing the region of the 
UV bump. Typical exposure times of about 10 hr led to an av- 
erage signal-to-noise ratio S/N of about 13 for the programme 
galaxies (see Noll et al. l2004t . The quality of the spectra is suf- 
ficiently high for the purpose of this study. The spectra used 
are S/N optimised co-added, flux calibrated (including a cor- 
rection for slit losses), and corrected for Galactic extinction 
{Eb-v = 0.018, Schlegel et al. using flie lTTMl MW 

extinction law. Finally, the atmospheric B and A absorption 
bands were removed from the observed spectra using an em- 
pirical band model. The final spectra are slightly contaminated 
by a second order spectrum for A > 6500 A. This effect does 
not significantly influence our study of the rest-frame UV con- 
tinuum; the resulting relative flux uncertainties are not greater 
than 2 to 3% in the UV-bump region for the selected galaxies. 

2.2. The lUE sample of nearby starbursts 

In order to test our new parametric description of the rest-frame 
UV SED of star-forming galaxies, we select a sample of 24 
well-known local starburst galaxies (z < 0.02; see Table |2j 
with available UV spectra (see Kinney et al. '1993' Calzetti 
et al. 1994 Heckman et al. 1998) . For this sample of nearby 
starbursts the luminosities at 1500 A are typically one order of 
magnitude fainter than those of the FDF galaxies (cf. Leitherer 
et al. .2002. and Noll et al. .2004j . 
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Table 2. The comparison sample of local starburst galaxies 
having lUE spectra with suitable S/N and spectral coverage. 
The (activity) class parameter shown is taken from Kinney 
et al. ST993i with the exception of IRAS 08339+6517 and 
NGC 1741, for which this parameter comes from Leitherer et 
al. (^002 1 and Johnson et al. J1999> . respectively. The mean- 
ing of the different parameters is as follows: BC(D)G - blue 
compact (dwarf) galaxy; Sy 2 = Seyfert 2 galaxy; SB nuc. = 
galaxy experiencing a starburst in its nucleus; Hs - hot spot 
galaxy; Lin = LINER; H II = galaxy with a spectrum typical 
of a H II region. All redshifts are taken from the NAS A/IPAC 
Extragalactic Database (NED). Values of Eb-v were calculated 
using the maps of Galactic dust IR emission of Schlegel et al. 
( fT998t . The oxygen abundances 0/H = 12 + log(0/H) are 
taken from Heckman et al. ( 1998 1 and complemented by data 
from Leitherer et al. (I2002> . Calzetti et al. (1994.) . and Johnson 
et al. ( fT^ . 



Name 


Class 


z 


Eb-v 
[mag] 


0/H 


ESO 338-4 


BCG 


0.0096 


0.09 


8.1 


IZwl8 


BCDG 


0.0025 


0.03 


7.2 


IRAS 08339-H6517 


HII 


0.0191 


0.09 


8.7 


Mrk33 


BCDG 


0.0049 


0.01 


8.4 


NGC 1097 


Hs + Lin 


0.0043 


0.03 


9.3 


NGC 1313 


HII 


0.0016 


0.11 


8.4 


NGC 1510 


BCDG 


0.0030 


0.01 


7.9 


NGC 1705 


BCDG 


0.0021 


0.01 


8.0 


NGC 1741 


nn 


0.0137 


0.05 


8.1 


NGC 2782 


SB nuc. 


0.0085 


0.02 


8.8 


NGC 3049 


SB nuc. 


0.0050 


0.04 


9.1 


NGC 3353 


BCDG 


0.0031 


0.01 


8.4 


NGC 3738 


HII 


0.0008 


0.01 


8.4 


NGC 42 14 


SB nuc. 


0.0010 


0.02 


8.2 


NGC 4385 


SB nuc. 


0.0071 


0.02 


8.7 


NGC 4449 


HII 


0.0007 


0.02 


8.4 


NGC 4670 


BCDG 


0.0036 


0.01 


8.2 


NGC 4861 


BCDG 


0.0028 


0.01 


8.0 


NGC 6052 


SB nuc. 


0.0157 


0.08 


8.6 


NGC 7496 


Sy 2 + H II 


0.0055 


0.01 


9.0 


NGC 7673 


HII 


0.0114 


0.04 


8.5 


NGC 7714 


SB nuc. 


0.0093 


0.05 


8.7 


NGC 7793 


HII 


0.0008 


0.02 


8.6 


UGC 9560 


BCDG 


0.0041 


0.01 


8.2 



The UV spectra of this reference spectroscopic sample are 
taken from the INES archive of lUE newly-extracted spectra. 
The resolution element of the low-resolution lUE spectra is 
about 6 A, i.e. similar to that achieved for the rest-frame UV 
spectra of the FDF sample of z ~ 2 UV-luminous galaxies 
(see Sect. 12. 1> . A fine adjustment of the resolution of the lUE 
spectra was not needed, since this study focuses on the analy- 
sis of the UV continuum. Further selection criteria are as fol- 
lows. First, all galaxies with a dominant young stellar popula- 
tion and available large 10" x 20" aperture, low-resolution UV 
spectra at short and long wavelengths were selected. In fact, 
only the combined spectrum obtained from the spectra taken 
with the short-wavelength prime camera (SWP) and the long- 



Table 3. Properties of the composite spectra of the FDF 2 < 
z < 2.5 and the lUE local galaxy sample in comparison. 
The table lists the number of spectra averaged, the mean red- 
shift, the mean equivalent width of six low-ionisation inter- 
stellar absorption lines (see Noll et al. I2004> . and the EWs 
of C IV AA 1548, 1550, He II A 1640, and C III] AA 1907, 1909. 
Emission lines have negative EWs. All errors given are mean 
errors. 



Sample 


FDF 


lUE 


A' 


34 


24 


z 


2.32 ± 0.02 


0.006 ± 0.001 


Wus [A] 


1.97 ±0.09 


1.89 ±0.08 


H^civ [A] 


3.80 ±0.16 


3.46 ± 0.37 


Wneii [A] 


-1.42 ±0.16 


-2.06 ± 0.45 


Weill] [A] 


-1.09 ±0.18 


-1.91 ±0.23 



wavelength prime and redundant cameras (LWP and/or LWR) 
cover the wavelength range between 1 150 and 3350 A, bracket- 
ing the region of the 2175 A feature. Second, only those galax- 
ies with spectra of suitable quality all over the UV spectral 
range were selected further. For obvious reasons, the data qual- 
ity at the overlapping region of the two lUE spectral channels 
at around 2000 A turned out to be a crucial parameter Hence, 
combined spectra exhibiting a break in this overlapping region, 
usually with rather poor S/N, were rejected. As a result, the 
combined spectra of the 24 selected lUE starbursts exhibit an 
average S/N of about 10 between 2000 and 2350 A, the final 
spectrum of each object being composed of all single spec- 
tra of sufficient S/N. In particular, the flux levels in the short 
and long wavelength spectra were adjusted in the wavelength 
range between 1850 A and 1980 A. Finally, the combined spec- 
tra were corrected for Galactic extinction and mapped into the 
rest frame. 

The correction for Galactic extinction makes use of the 
maps of Galactic dust emission at IR wavelengths (Schlegel 
et al.'1998^ for the derivation of realistic values for Eb-v, and 
of the CCM Milky Way extinction law for Ry - 3.1. This step 
is crucial, since the IVIW extinction law for the average ISM ex- 
hibits a relatively strong UV bump, so that small errors in Eb-v 
can lead to important uncertainties in the evaluation of the UV 
bump strength from the spectra of the 24 selected lUE star- 
bursts. However, for the sample galaxies, Eb-v ranges between 
0.01 and 0.11, with an average value of 0.03 (see Table|2j, so 
that the effect of the correction for Galactic extinction on most 
of the spectra is possibly minimised. 

2.3. Basic properties of the FDF and lUE samples 

Fig. 121 reproduces the composite spectrum of the 34 UV- 
luminous FDF galaxies at 2 < z < 2.5 and the composite 
spectrum of the 24 local lUE starbursts. These composite spec- 
tra were calculated according to the recipe given in Noll et 
al. J2004t . This means that the rest-frame spectra were first 
normalised to a common mode in the range 1250 to 1500 A 
and then averaged using equal weights. Wavelength regions af- 
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o 

CO 

oi 



1500 2000 
Wavelength [A] 

Fig. 2. Composite of 34 FDF galaxy spectra with S/N > 4 in the 
range 2 < z < 2.5 (upper spectrum in the upper panel) in com- 
parison with a composite of 24 lUE spectra of local galaxies 
(lower spectrum in the upper panel; offset by -0.5 for clarity). 
The lower panel gives the ratio of the FDF and the lUE com- 
posite spectrum. The positions of prominent UV lines are indi- 
cated. The Lya line of the lUE composite is strongly affected 
by emission from the geocorona. 



fected by residuals of strong night sky lines or other defects 
were excluded from the averaging. 

The lower panel of Fig. |2 shows the ratio of the FDF 
and lUE composites. In general, both rest-frame UV spectra 
are very similar, and show clear signatures of the dominance 
of young stellar populations. For a large wavelength range 
the continuum flux increases towards shorter wavelengths. 
However, for A < 1800 A the FDF composite is distinctly flat- 
ter than the lUE composite. Conversely, for A > 1800 A the 
local starburst composite shows a flatter slope than the other 
We will discuss these behaviours in the following sections. 

In addition to the overall UV continua, the average equiv- 
alent width of six prominent low-ionisation lines (Wlis, see 
Shapley et al. 2003 Noll et al.p004) are similar between the 
two samples (Fig.Hand Table|3ll. In fact, Wus = 1 .97 ± 0.09 A 
for the FDF composite and 1.89 + 0.08 A for the local one. 
The quoted errors include statistical uncertainties and the vari- 
ance of the line strengths within each sample (see Noll et al. 
[2004 1. These low-ionisation lines originate in the intervening 
neutral interstellar medium between the OB stars (responsi- 
ble of the rest-frame UV continuum) and the observer. Due to 
their high saturation, these lines are barely sensitive to chemi- 
cal abundances. Hence, they primarily trace the spatial and ve- 
locity distribution of the neutral interstellar gas illuminated by 
the rest-frame UV continuum. Therefore, Wlis is a measure of 
the product of the neutral clouds' covering fraction and the ve- 
locity spread. According to Shapley et al. (2003 ) differences in 
Wlis are mainly due to variations in the covering fraction. 



Comparable equivalent widths are also obtained for the im- 
portant high-ionisation C IV doublet: 3.80 ± 0. 14 A for the FDF 
sample, and 3.46 + 0.37 A for the lUE sample. The C IV equiv- 
alent width mainly measures the strength of winds from hot, 
massive and, thus, young stars and is a good metallicity indi- 
cator (Leitherer et al. 2001 1. Using C IV Mehlert et al. HM^ ^ 
estimated an average metallicity of Z 0.42 for an analo- 
gous sample of FDF galaxies at 2 < z < 2.5. Measuring the 
purely photometric indices "1370" and "1425" (Leitherer et al. 
i2001j in medium-resolution spectra of two z ~ 2 FDF galaxies 
Mehlert et al. ^005 1 find indications for even higher metallic- 
ities close to the solar value. A similar metallicity as the C IV 
based of the FDF galaxies (assuming the solar oxygen abun- 
dance (0/H)q = 8.93) we find for our lUE sample of nearby 
starbursts having (O/H) w 8.6 (see Tableland, e.g., Heckman 
etal. 1998J . 

In addition to the P Cygni profiles of high-ionisation lines 
such as CIV, the detection of the Hell and CIII] emission 
features indicates the presence of a significant population of 
young, very massive stars in both samples (Fig.|2land Table|3}. 
In particular, the He II line emission is essentially produced by 
Wolf-Rayet stars (Schaerer i2003il . whereas the C III] line emis- 
sion has a nebular origin, and traces the electron temperature 
of H II regions (see Heckman et al. ll998> . 

3. Analysis 

In this study we compare diagnostics measured on the rest- 
frame UV spectra of the 34 FDF sample galaxies (see Sect. l2.1> 
with those determined on the rest-frame UV SEDs of models of 
starburst galaxies. This procedure is analogous to that adopted 
by Gordon et al. J1997> for the analysis of the dust proper- 
ties for the Calzetti et al. (1994'( sample of nearby starbursts. 
In addition, we use the lUE sample of 24 local starbursts (see 
Sect. 12. 2> as a further comparison. Here we illustrate both the 
models and the diagnostics adopted in this study. Furthermore, 
we demonstrate that the new diagnostics tailored to the FORS 
data lead to the same conclusions as those of Gordon et al. 



1997> when applied to the lUE sample. 



3.1. Modelling the UV continuum of starburst galaxies 

We produce synthetic SEDs for starburst galaxies by combin- 
ing stellar population synthesis models with models of radiative 
transfer of the stellar and scattered radiation through the dusty 
ISM. 



3.1 .1 . Stellar population synthesis models 

Intrinsic (i.e. non attenuated by dust) SEDs for starburst galax- 
ies are derived from the stellar population synthesis code for 
the stellar continuum emission of Maraston (2005, .MOSl ^. 

^ The EWs of C IV listed in Mehlert et al. <2002> are on average 
about 20% larger than in this study. This is due to a wider wavelength 
range to derive the continuum level (190 A instead of ~ 80 A) and an 
automatic instead of interactive measurement technique. 

^ The Maraston (2005) models are available at 
http : //www-astro . physics . ox . ac . uk/ ~maraston/. 
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These models include the thermally pulsating asymptotic gi- 
ant branch phase for intermediate-mass (2 - SM©) stars, at 
variance with all others (e.g. Fioc & Rocca-Volmerange 1997, 
IFR97 Leitherer et al. 1999 Bruzual & Chariot 2003. .BC03« . 
However, this relevant difference does not affect any of the 
results reported here, since the SEDs produced for analogous 
models of M05. .BC03i and iFR9 7 do not differ significantly in 
shape at rest-frame UV wavelengths (see, e.g.. Fig. 14 in M05 1. 
This is also shown in Table |5] of Sect. 13.2.21 Finally, we note 
that the spectral resolution of the M05. models (AA - 10 A) 
matches reasonably well those of the rest-frame UV continua 
of the FDF sample galaxies, as given by the FORS spectra 
(AA - 5 - 8 A), and of the nearby starbursts observed with 
lUE (A^ = 6 A). 

In Sect. 12.31 we have shown that there is evidence for the 
z ~ 2 UV-luminous FDF galaxies being both actively forming 
stars and metal enriched ((Z) > 0.4 Zq). In addition, starburst 
galaxies in the nearby and distant universe seem to undergo 
multiple star-formation episodes of a few lOMyr to several 
lOOMyr (Calzetti .1997. Shapley et al. 120011 van Dokkum et 
al. I2004> . Therefore, we adopt composite stellar populations 
models where the SFR is constant with time, the metallicity 
is fixed and ranges between 2% and twice solar, and the time 
elapsed since the start of the star-formation activity (i.e. the age 
of the model) ranges between 5 and 500 Myr As a reference, 
we compute additional simple stellar population models span- 
ning the same metallicity range as for the previous models but 
with ages between 5 and 100 Myr. All these models assume a 
Salpeter (.1955 J initial mass function with an upper mass cut-off 
of 100 Mo. 

Given the measured metallicities (Mehlert et al. 2002 
12005 't and the expected ages of the z ~ 2 UV-luminous galax- 
ies, in the following analysis we will make use mostly of M05l 
models with a continuous star-formation activity (i.e. with a 
SFR e-folding time r = lOGyr) along 50, 200, or 500 Myr, 
and metallicities equal to 0.5 or l.OZ©. These models are 
also suitable to reproduce the properties of the lUE sample of 
nearby starbursts, though a few objects in the latter sample have 
Z < 0.1 Zq (e.g. IZwl8) and a few others have significantly 
supersolar metallicity (see Table|2li. 

3.1 .2. Radiative transfer models 

We make use of the radiative transfer calculations for the 
SHELL dust/stars configuration of lWGOOl in order to describe 
the properties of dust attenuation for the galaxies under inves- 
tigation. The reason is that the rather regular properties of dust 
attenuation in starburst galaxies of the local universe (Calzetti 
et al. 119941 12000") may be reproduced by a model where a 
two-phase, clumpy, dusty medium is distributed in a shell sur- 
rounding the stars (Gordon et al. ll997> . This dust/stars config- 
uration is analogous to that adopted by Fischera et al. (2003) 
to give an independent physical explanation of the so-called 
"Calzetti law" (Calzetti et al. [T994ll2000t based on a turbulent 
ISM. According to Dopita et al. (2005), the pressure in the dif- 
fuse phase of the ISM (or, equivalently, its density), and the 
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Fig. 3. Ratio of the attenuation factor A at the wavelength A to 
the average attenuation in the V band as function of A for dif- 
ferent extinction laws. The ICCMl Milkv Way extinction law is 
shown for Ry - AvIEb-v = 3.1 as solid line. The LMC aver- 
age, LMC 2 supershell, and SMC bar extinction laws (Gordon 
et al. '2003 1 are marked by dash-dotted, dashed, and dotted 
lines, respectively. 

molecular-cloud dissipation timescale control the pan-spectral 
SED of starburst galaxies. 

The tWGOO radiative transfer models (including multiple 
scattering) are based on Monte Carlo calculations for a scale- 
free spherically-symmetric geometry (see also Gordon et al. 
1200 U . In particular, the SHELL dust/stars configuration as- 
sumes the stars to be distributed inside an inner sphere with 
radius equal to 0.3 times the system radius, and the dust to 
be distributed inside a shell extending radially from 0.3 to 1 
times the system radius. Dust can be distributed locally either 
in a homogeneous (i.e. diffuse) ISM, or in a two-phase, clumpy 
ISM. In the latter case, high-density clumps are embedded in a 
low-density interclump medium, the low-to-high-density ratio 
being equal to 0.01, with a filling factor of 0.15. Models are 
parameterised by the radial extinction optical depth from the 
centre to the edge of the dust environment at V band (ry), as- 
suming a constant density homogeneous distribution, tv rang- 
ing from 0.25 to 50. Finally, the mixture of dust grains can be 
characterised by an extinction curve that is typical either of the 
SMC or of the MW (see Fig.|3land Sect.Q. 

Hereafter, we illustrate a few aspects of these SHELL ra- 
diative transfer models for the non-specialist reader. Fig. |4] 
shows the behaviour of the UV spectral slope as defined 
in Sect. 13.2.11 as a function of the attenuation (in mag) at 
1500 A A 1500 for different starburst models with age equal to 

* The non-ionising UV-flux density distribution of an observed star- 
forming galaxy is well reproduced by a power law ~ X^^, where /3 is 
referred to as the UV slope (see, e.g., Calzetti et al. ll994t . 
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Fig. 4. The UV continuum slope as a function of the at- 
tenuation factor at 1500 A A 1500 for different reddened model 
spectra. In all cases a M05 model with Salpeter IMF, contin- 
uous star-formation, an age of 200 Myr, and solar metallicity 
was combined with SHELL models of lWGOOl for different dust 
distribution, extinction laws, and optical depth. Different dust 
column densities are indicated by filled squares (MW dust) 
and lozenges (SMC dust), connected with solid (homogeneous 
model) and dashed lines (clumpy model), respectively. The 
symbols have been plotted for the average radial visual opti- 
cal depth TV = 0, 0.25, 1, L5, 6, 7, 10. The dash-dotted 
line refers to the Calzetti et al. 



( I2000t attenuation law. The 
crosses are plotted in intervals of AEb-v = 0.1. 



200 Myr and solar metallicity (cf. Sect. 13. LH . The relation be- 
tween jS and A 1500 is linear in case of a homogeneous dusty 
ISM, whatever the dust type. Conversely, it is non-linear for 
the two-phase, clumpy dust distribution. In the former case, /3 
and A 1500 increase at a constant rate because the angle-averaged 
optical depth at UV wavelengths increases monotonically with 
the dust column density (i.e. Ty). In the latter case, first /? and 
A 1500 show a positive coiTelation because the optical depth at 
UV wavelengths of a clump increases with Ty- However, when 
the clumps become optically thick at UV wavelengths, the at- 
tenuation at these wavelengths is given by the simple block- 
ing action of the clumps, while the UV radiation reaches the 
outside observer through the optically thin, diffuse component 
of the ISM. Hence, /3 stalls while A 1500 keeps on increasing. 
When the angle-averaged optical depth at UV wavelengths of 
the diffuse component of the ISM tends to unity, and A 1500 
increase together again. We note that, in case of a more real- 
istic size and density distribution of the dusty clumps, one can 
expect a behaviour that is intermediate between those produced 
by the homogeneous and two-phase, clumpy dust distributions. 
Furthermore, Fig. |3 confirms that SMC-type dust produces a 



Table 4. Wavelength windows for the derivation of the UV con- 
tinuum slope parameters /3 and 71 to 74 using power-law fits 

Par. Windows 



/3 


1268- 


1284, 1309 - 


1316, 1342 - 


1371, 1435 - 


1496, 




1562- 


1583, 1677 - 


1740 






71 


1268 - 


1284, 1309 - 


1316, 1342 - 


1371, 1435 - 


1496 


72 


1562- 


1583, 1677 - 


1740, 1760 - 


1833, 1866- 


1890 


73 


1920- 


2050, 2070 - 


2175 






74 


2175 - 


2310, 2400- 


2480 
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Fig. 5. Demonstration of the derivation of the UV continuum 
slopes /3 and 71 to 74 for the composite spectrum of 34 lumi- 
nous FDF galaxies at 2 < z < 2.5 (cf. Fig.|2}. The power law 
fits are shown as thick dotted (J3) and dashed (71 to 74) curves. 
The limits for the fit ranges are indicated by vertical dotted (J3) 
and dashed (71 to 74) Unes. 



larger attenuation at UV wavelengths than MW-type dust, for 
fixed Ty, whatever the structure of the dusty ISM. 

Finally, Fig. |4] shows that a linear relation between y6 and 
A 1500 is expected when the empirical "Calzetti law" is assumed 
to hold, whatever the dust column density of the ISM. We re- 
mind that WGOOi could achieve the best agreement with the 
"Calzetti law" for a two-phase, clumpy SHELL model with 
SMC dust, and ry = 1.5. 

3.2. A new parametric description of the UV continuum 
3.2.1 . The parameter set 

The UV SEDs of composite stellar population models can be 
parameterised with a power law (/(/I) oc A^), since the bulk ra- 
diation comes from the emission from hot stars in the Rayleigh- 
Jeans regime. This is so also true for observed UV SEDs of 
nearby starbursts (e.g. Calzetti et al. .l994j . that are affected by 
dust attenuation (see footnote 3). The latter authors measured 
the UV slope yS across the wavelength range 1200 - 2600 A, 
with the exclusion of the wavelength region between 1950 and 
2400 A. By doing so, their /? was not affected by the possible 
presence of dust characterised by an extinction law with the 
UV bump. Other analyses of the UV continuum properties of 
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star-forming galaxies in the local and distant universe defined 
P across a narrower range, i.e. from 1200 to 1800 A (Heckman 
et al l 19981 Leitherer et al. 2002 Noll et al. 2004 1, mostly ow- 
ing to the missing spectral coverage at rest-frame A > 2400 A. 
Following the previous authors we also measure j3 across the 
wavelength range 1200 - 1800 A. The quaUty of the FORS 
spectra available for the 34 FDF galaxies at 2 < z < 2.5 is 
not good enough to provide a reliable fit at /I > 2400 A. 

In addition to the "standard" parameter/?, we define six new 
parameters, since the main interest of this study is to investi- 
gate if and when the UV bump is present among the FDF sam- 
ple galaxies. The definition of four parameters (see Table|3i is 
based on power-law fits to sub-regions of the UV continuum as 
well as that of In particular, we divide the wavelength range 
between 1200 and 2500 A in four sub-regions, delimited by the 
CIV doublet at 1550 A, the CIII] emission feature at 1908 A, 
and the centre of the UV bump at 2175 A. Furthermore, we se- 
lect windows within a given sub-region (see Table 13 in order 
to avoid that the presence of strong spectral lines produces a 
bias of the fit to the continuum. For A < 1800 A we essen- 
tially adopt the windows defined by Calzetti et al. (Q994 1, with 
the exception of the 1435 - 1496 A window, which is taken 
from Leitherer et al. ^)Q2\. We call these four parameters yi, 
72, 73, and 74 after the progressive number of the four sub- 
regions, starting from that defined by the shortest wavelengths. 
Fig. |5]reproduces the five power-law fits corresponding to the 
previous five parameters (i.e. yS plus the four new slopes of 
the fits) for the composite spectrum of 34 UV-luminous FDF 
galaxies at z ~ 2. In addition, we introduce two parameters, 
i-e. 7i2 =71-72 and 734 = 73 - 74. The meaning of these 
two additional parameters is straightforward: 712 measures the 
"curvature" of the spectral UV continuum across 1500 A; 734 
measures the "curvature" of the spectral UV continuum across 
the UV bump (negative 734 identify the presence of dust with 
an extinction law exhibiting the 2175 A feature). The parameter 
734 does not represent the "strength" of the UV bump however. 
In fact, radiative transfer effects tend to dilute the prominence 
of the UV bump in the resulting attenuation function with re- 
spect to that in the extinction curve when the opacity increases, 
whatever the dust/stars configuration and the structure of the 
dusty ISM, with the exception of a homogeneous SHELL con- 
figuration (Cimatti et aL nWHIWHUHl Pierini et al. EUn^t . 

3.2.2. Sensitivity of tine parameters to properties of tine 
stellar populations or of dust attenuation 

The five parameters /3, and 71 to 74 describe the observed UV 
SEDs of a galaxy and, thus, depend in a complex way on the 
properties of the stellar populations and of dust attenuation (see 
Sect. 13. 1> . Nevertheless, we can disentangle the principal de- 
pendencies of each parameter by investigating the behaviour 
of the five slopes, as obtained for the UV SEDs of a suite of 
starburst models, as a function of the model parameters. 

In Fig.|6lwe illustrate the effects on the UV continuum pa- 
rameters /3 and 73 of diff'erent metallicities and ages, for single- 
burst and continuous star-formation models of 'M05" A pri- 
ori, one would expect that 73 is particularly sensitive to very 
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Fig. 6. The relation between the slope parameters 73 and yS 
for diff'erent IM05I model spectra. The tracks plotted show the 
change of /3 and 73 with increasing age (f < 1 Gyr) for stellar 
population models diff'ering in the star-formation history (sin- 
gle burst or continuous star-formation) as well as the metallic- 
ity (0.02, 0.5, 1.0, and 2.0Zo). For single burst models (dotted 
lines) symbols (see legend) are plotted for 5, 20, and 100 Myr. 
Ages of 5, 50, and 500 Myr are marked in case of models 
with continuous star-formation (solid lines). All tracks start at 
y6 3 and 73 3. 

young ages and high metallicities, since it characterises the 
wavelength range from 1900 to 2175 A. However, 73 hardly 
changes for reasonable models of z ~ 2 UV-luminous galax- 
ies, with continuous star-formation, metallicity equal to 0.5 or 
I.OZ0, and ages higher than a few lOMyr (see Sect. B.l.Tt . The 
same conclusion holds for /3: conspicuous changes in /3 would 
be possible only for single burst models with ages greater than 
20 Myr The latter models are not realistic, since they do not 
provide large UV luminosities long after the episode of star- 
formation activity. Furthermore, values of > -2 are produced 
by stellar population synthesis models with a SFR e-folding 
time as short as 1 Gyr only about 10 Gyr after the onset of the 
star-formation activity, i.e. at a time much longer than the age 
of the universe at z ~ 2 (i.e. ~ 3 Gyr). SFR e-folding times less 
than 1 Gyr can not justify the vigorous star-formation activity 
of the UV-luminous FDF galaxies at z ~ 2 that is witnessed by 
the measured emission Une EWs (see Tabled}- 

Table|5]lists typical values (see Sect. B.l.Tt for all five UV 
continuum slope parameters. In particular, it shows that 74 ex- 
hibits only a minor dependence on the properties of the stel- 
lar populations, while 72 is particularly sensitive to changes in 
metallicity and age. The latter is explained by the great num- 
ber of spectral lines originated in the atmospheres of hot stars 
that fall within the wavelength range 1550 to 1900 A, where 72 
is defined, and the sensitivity of the same wavelength range to 
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Table 5. Continuum slope parameters for continuous star- 
formation models of M05 differing in metallicity (0.5 and 
1 Zq) and age (50 and 500 Myr), and for analogous models of 
IBC03l and .FR97. 



Par. 


0.5 Zo 


IZo 


0.5 Zo 


IZo 


Ref. 


P 


-2.74 


-2.69 


-2.49 


-2.44 


M05 


7i 


Z.oU 


1 9.1 

— Z.o / 








72 


-2.23 


-1.98 


-2.00 


-1.74 


M05 


73 


-2.16 


-1.93 


-2.04 


-1.85 


M05 


74 


-3.04 


-3.03 


-2.92 


-2.96 


M05 


p 


-2.73 


-2.67 


-2.45 


-2.39 


BC03 


7i 


-2.83 


-2.96 


-2.40 


-2.50 


BC03 


72 


-2.23 


-1.86 


-1.97 


-1.60 


BC03 


73 


-2.07 


-1.75 


-1.95 


-1.68 


BC03 


74 


-3.00 


-3.01 


-2.88 


-2.95 


BC03 


p 


-2.61 


-2.58 


-2.36 


-2.33 


FR97 


7i 


-2.58 


-2.70 


-2.19 


-2.31 


FR97 


72 


-2.24 


-1.94 


-2.01 


-1.71 


FR97 


73 


-2.11 


-1.85 


-2.00 


-1.79 


FR97 


74 


-2.85 


-2.88 


-2.76 


-2.84 


FR97 
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Fig. 7. The relation between the slope parameters 73 and ji for 
different reddened model spectra. The hne types and symbols 
are as in Fig.|3 

the O-to-B-star number ratio. Finally, Table |5] shows that the 
slope parameters span a range Ay of up to 0.5 in response to 
the variations in the properties of the stellar populations that 
are assumed to be reasonable for our sample galaxies. On the 
other hand, Table|5]shows that the different y parameters char- 
acterising the rest-frame UV SED can change by up to 0.2 with 
respect to those obtained from the IM05I models, when they 
are determined from intrinsic rest-frame UV SEDs taken from 
[BC03i Or .FR97. Hence, the dependence of all the seven UV con- 



tinuum slope parameters on the source of the stellar population 
synthesis code does not affect the study of the extinction law in 
UV luminous galaxies at z ~ 2. 

Now, stellar population synthesis models do not include the 
absorption of UV light by interstellar gas, so that one may 
expect that the parameters listed in Table |5l are affected by 
systematic uncertainties. Therefore, we have computed semi- 
analytic spectra including line absorption by the ISM by com- 
bining the normalised, composite rest-frame UV spectrum of 
our high-redshift sample galaxies (Fig.O with the UV spectra 
of the M05 models listed in Table |5l As a result, we estimate 
that the deviations in the slope parameters listed there, as pro- 
duced by interstellar line-absorption, are typically about 0.1. 
This very small effect is also a consequence of the careful selec- 
tion of wavelength windows excluding strong interstellar lines, 
where the five UV continuum slope parameters are defined (see 

Fig. 13. 

Finally, we discuss how the latter parameters depend on 
the properties of dust attenuation, for the models discussed 
in Sect. 13.1.21 In particular, we illustrate the behaviour of 
dusty, star-forming galaxies with a solar metallicity, 200 Myr- 
old composite stellar population in the y6 - 73 plane (see Fig.0 
already introduced (see Fig. |6}. In this diagram, the loci of 
models with MW and SMC dust are clearly separated, what- 
ever Ty. The slope parameters 73 and fi of the two models differ 
by more than 1.6 and 0.5, respectively, already for ^^^ > 0.25. 

Furthermore, 73 2 is a delimiter of the regions populated by 

the models for MW-like and SMC-like dust. Hence, the distri- 
bution in the J6-73 plane can reveal the presence of dust with an 
extinction law exhibiting the UV bump already for low values 
of the dust column density. On the other hand, it is not particu- 
larly sensitive to the structure of the dusty ISM. As a reference, 
Fig.0reproduces additional models with the same stellar pop- 
ulations as the others but where dust attenuation follows the 
Calzetti et al. ( 2000 1 attenuation law. As expected (Gordon et 
al. II997> . the latter models show a similar behaviour to those 
with a SHELL configuration of SMC-type dust distributed in a 
two-phase, clumpy medium. 

In conclusion, the five UV continuum slopes defined in this 
study are proved to be much more sensitive to the character- 
istic extinction law (i.e. to the properties of the mixture of the 
dust grains) of the system, than to those of the stellar popula- 
tions present in the system, for a fixed dust/stars distribution. 
Conversely, they cannot constrain the structure of the dusty 
ISM. 

3.2.3. A test for the method: the IDE sample of nearby 
starbursts 

Here we test the reliability of our method of investigating the 
presence of dust with an extinction law exhibiting the 2175 A 
bump, that is based on a new parametric approach. We do so 
by investigating the slope parameters defined in Sect. l3.2.T1 of 
our lUE sample of nearby starbursts (see Sect. 12. 2> with re- 
spect to models. Table |6l lists the parameters measured by us 
for this reference sample. Composite stellar populations mod- 
els are calculated as in .M05i for a standard Salpeter IMF, a 
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Table 6. UV continuum slope parameters of our lUE sample of 24 local starburst galaxies (see Table|2J 



Name 


P 


71 


72 


73 


74 


712 


734 



ESO 338-4 


-1, 


,19 


± 


0, 


,09 


+0, 


,02 


+ 


0, 


.08 


-1. 


.06 


+ 


0, 


.07 


-1 


.55: 


tO. 


.09 


-3. 


.01 


± 


0, 


.07 


+ 1. 


,08 


+ 


0, 


.11 


+ 1, 


.46: 


tO. 


.11 


IZwl8 


-1. 


,80 


± 


0, 


,07 


-1, 


,20 


+ 





,06 


-2. 


.05 


+ 


0, 


.06 


-2, 


.91 : 


tO. 


.23 


-3. 


.83 


+ 


0, 


.22 


+0. 


,86 


+ 





.08 


+0, 


.92 ± 0, 


.32 


IRAS 08339+6517 


-0. 


,63 


± 


0, 


,08 


+0, 


,16 


± 





,07 


-1. 


.45 


± 


0, 


.09 


-2, 


.15: 


tO. 


.14 


+0. 


.29 


+ 


0, 


.13 


+ 1. 


,61 


+ 





.11 


-2, 


.44 ± 0. 


.19 


Mrk33 


-0. 


,95 


± 


0, 


,09 


-1, 


,30 


± 


0, 


.07 


-0. 


.48 


+ 





.11 


-1, 


.39: 


tO. 


.14 


-0. 


.76 


+ 





.07 


-0. 


,82 


+ 





.13 


-0 


.63 : 


t 0. 


.16 


NGC 1097 


-0. 


,37 


± 


0, 


,13 


-0, 


,75 


+ 





,12 


+0. 


.77 


+ 





.12 


+0, 


.91 : 


tO. 


.17 


-0. 


.65 


+ 


0, 


.15 


-1. 


,51 


+ 





.17 


+ 1, 


.56 : 


t 0. 


.23 


NGC1313 


-1. 


,21 


± 


0, 


,19 


-0, 


,58 


+ 





.20 


-0. 


.87 


+ 





.14 


-3, 


.33: 


tO. 


.28 


-1. 


.72 


+ 


0, 


.22 


+0. 


,29 


+ 





.24 


-1, 


.61 : 


t 0. 


.36 


NGC 1510 


-1. 


,25 


+ 


0, 


,08 


-0, 


,62 


+ 





.07 


-0. 


.92 


± 


0, 


.08 


-1, 


.54 ± 0, 


.17 


-1, 


.74 


+ 


0, 


.14 


+0. 


,30 


± 


0, 


.10 


+0, 


.21 : 


t 0. 


.22 


NGC 1705 


-2. 


36 


+ 


0, 


,04 


-2, 


,10 


+ 


0, 


.04 


-2. 


.13 


+ 


0, 


.06 


-1 


.91 ± 0, 


.05 


-3. 


,56 


+ 


0, 


.05 


+0. 


,03 


+ 


0, 


.07 


+1 


.65 : 


tO. 


.07 


NGC 1741 


-0. 


,16 


+ 


0, 


,13 


+ 1, 


,15 


+ 


0, 


.11 


-0, 


.64 


+ 


0, 


.11 


-0, 


.71 : 


tO. 


.14 


-2. 


,58 


+ 


0, 


.18 


+1. 


,79 


+ 


0, 


.16 


+1, 


.87 : 


t 0. 


.23 


NGC 2782 


-0. 


,79 


+ 


0, 


,12 


-0, 


,36 


+ 


0, 


.11 


-0. 


.50 


+ 


0, 


.10 


+0, 


.31 : 


tO. 


.24 


-0. 


,07 


+ 


0, 


.26 


+0, 


,14 


+ 


0, 


.14 


+0, 


.38 : 


tO. 


.35 


NGC 3049 


-1. 


,23 


+ 


0, 


,12 


-0, 


,65 


+ 


0, 


.08 


-0. 


.43 


+ 


0, 


.12 


-2, 


.20: 


tO. 


.27 


-0. 


,94 


+ 


0, 


.17 


-0. 


,22 


+ 


0, 


.14 


-1 


.27 : 


tO. 


.32 


NGC 3353 


-0, 


,84 


± 


0, 


,07 


-0, 


,37 


± 


0, 


.07 


-1. 


.42 


± 


0, 


.06 


-2, 


.50: 


tO. 


.21 


-1. 


,68 


+ 


0, 


.15 


+1. 


,05 


+ 


0, 


.09 


-0, 


.82 : 


tO. 


.26 


NGC 3738 


-1, 


,47 


± 


0, 


,07 


-0, 


,94 


+ 





.06 


-1. 


.11 


± 


0, 


.05 


-0, 


.97: 


tO. 


.14 


-2. 


,06 


+ 


0, 


.12 


+0. 


,17 


+ 


0, 


.08 


+1, 


.09: 


to. 


.18 


NGC 4214 


-0, 


,98 


± 


0, 


,05 


-0, 


,24 


+ 


0, 


.04 


-1. 


.20 


+ 


0, 


.03 


-1 


.07: 


tO. 


.07 


-2. 


,35 


+ 


0, 


.06 


+0. 


,97 


+ 


0, 


.05 


+1, 


.29: 


to. 


.10 


NGC 4385 


-0. 


,16 


± 


0, 


,12 


+0, 


,10 


+ 


0, 


.12 


-0. 


.88 


+ 


0, 


.09 


+0, 


.06: 


tO. 


.14 


-2. 


,88 


+ 


0, 


.12 


+0, 


,98 


+ 


0, 


.15 


+2, 


.94: 


to. 


.18 


NGC 4449 


-1, 


,31 


± 


0, 


,05 


-0, 


,56 


± 


0, 


.03 


-1. 


.19 


+ 


0, 


.04 


-1 


.21 : 


tO. 


.09 


-3. 


,56 


+ 


0, 


.09 


+0, 


,63 


+ 


0, 


.05 


+2, 


.35: 


to. 


.12 


NGC 4670 


-1. 


,23 


+ 


0, 


,06 


-0, 


,97 


+ 


0, 


.05 


-1. 


.45 


+ 


0, 


.05 


-0, 


.21 : 


tO, 


.09 


-1. 


,77 


± 


0, 


.07 


+0, 


,49 


+ 


0, 


.07 


+ 1 


.56: 


to. 


.12 


NGC 4861 


-2, 


,03 


± 


0, 


,05 


-1, 


,54 


+ 


0, 


.05 


-2. 


.00 


+ 


0, 


.05 


-2, 


.28: 


tO. 


.07 


-2. 


,30 


+ 


0, 


.07 


+0, 


,45 


+ 


0, 


.07 


+0, 


.01 : 


to. 


.10 


NGC 6052 


+0. 


,39 


± 


0, 


,18 


+2, 


,38 


+ 


0, 


.15 


-0. 


.91 


+ 


0, 


.13 


-0, 


.53: 


tO. 


.18 


-1. 


,47 


+ 


0, 


.14 


+3, 


,29 


+ 


0, 


.20 


+0, 


.94: 


to. 


.23 


NGC 7496 


-1. 


,07 


+ 


0, 


,12 


-1, 


,26 


+ 


0, 


.10 


+0. 


.19 


+ 


0, 


.11 


-3, 


.71 : 


tO. 


.15 


-0. 


,97 


+ 


0, 


.10 


-1, 


,45 


+ 


0, 


.14 


-2, 


.74: 


to. 


.18 


NGC 7673 


-0. 


,96 


+ 


0, 


,09 


-0, 


,50 


± 


0, 


.09 


-0. 


.63 


+ 


0, 


.08 


-2, 


.29: 


tO. 


.16 


-2. 


,05 


+ 


0, 


.11 


+0, 


,13 


+ 


0, 


.12 


-0, 


.24: 


to. 


.20 


NGC 7714 


-0. 


,95 


+ 


0, 


,05 


-0, 


,90 


± 


0, 


.04 


-0. 


.37 


+ 


0, 


.05 


-2, 


.29: 


tO, 


.16 


-1. 


,24 


± 


0, 


.14 


-0, 


,52 


+ 


0, 


.06 


-1 


.05: 


to. 


.21 


NGC 7793 


-0. 


,69 


± 


0, 


,13 


+0, 


,43 


+ 





.13 


-1. 


.17 


+ 


0, 


.10 


-4, 


.42: 


tO. 


.19 


-2. 


,13 


+ 


0, 


.18 


+1, 


,60 


+ 


0, 


.16 


-2, 


.29: 


to. 


.26 


UGC9560 


-1, 


,99 


± 


0, 


,06 


-1, 


,74 


+ 





.05 


-1. 


.54 


± 


0, 


.07 


-1 


.27 ± 0, 


.12 


-2. 


,51 


+ 


0, 


.13 


-0, 


,20 


+ 


0, 


.08 


+ 1, 


.24: 


to. 


.18 



continuous star-formation activity along 50 or 500 Myr, and 
a metallicity equal to 0.5 Z© or 1 .0 Z©. Furthermore, the prop- 
erties of dust attenuation are described through the two-phase, 
clumpy SHELL dust/stars configuration of WGOO for SMC- 
like or MW-like dust, where the average opacity of the dust 
shell is a free parameter. 



The combined slope parameter 734 is most suitable to study 
the presence of the 2175 A feature. This is illustrated by the plot 
of versus y6 in Fig. |8l which shows that models with differ- 
ent composite stellar populations and SMC-type dust, lacking 
the 2175 A feature, exhibit 734 ~ 1, whatever Ty and respec- 
tively. Conversely, 734 rapidly decreases below zero when fi in- 
creases for models where UV-bump producing MW-type dust 
is present (see Fig.|3}- Now, the 24 local starbursts of the lUE 
sample are distributed in a region delimited by -2.5 < P < 0.5 
and -3 < 734 < 3. The bulk of these galaxies populates a re- 
gion of the p - 734 plane that is confined by the model curves, 
for yS > -3. Our interpretation is that these nearby starbursts 
possibly contain relatively small amount of dust with an extinc- 
tion law exhibiting only a modest UV bump, if any. Moreover, 
analysing the distributions of the other slope parameters (see 
Fig. |5) we find that the far-UV uprise of the extinction law is 
intermediate between those of the SMC and MW. These are 
the same conclusions of Gordon et al. ( 1997 1 and WG OO for 
the sample of 30 nearby starbursts with lUE data selected by 
Calzetti et al. J1994> . Hence, our method of analysis is proved 
to be reliable. 



4. Results 

Table lists the values of the UV continuum slopes yS and 71 
to 74, together with the parameters 712 and 734 measured for 
the FDF sample of 34 UV-luminous galaxies at 2 < z < 2.5 
(see Sect. 12. H . Hereafter, we will investigate the distribution 
of these galaxies in diagnostic diagrams defined by p and one 
of the six 7 parameters in comparison with that of the 24 local 
starbursts of the lUE sample (see Sect. l2.2l and Table|6} and of 
models combining stellar population synthesis and dust atten- 
uation as in Sect. 13.2.^ 

Meurer et al. ( 1999 1 proposed the parameter as a measure 
of dust attenuation at UV wavelengths for nearby and distant 
starburst galaxies. More recently, Leitherer et al. ( I2002> deter- 
mine 

Ai5oo[mag] = 2.19 (yS -;8o), = -2.5 (1) 

as the relation between p and the attenuation at 1500 A for 
nearby starburst galaxies. However, radiative transfer mod- 
elling ( IWGOO> demonstrates that the "standard" yS parameter is 
a model-dependent measure of dust attenuation at UV wave- 
lengths. In particular. Fig. 0] illustrates the relation between 
P and A 1500 for the WGOO models discussed in Sect. 13.1.21 
Furthermore, it shows that yS is not a linear measure of dust 
attenuation at UV wavelengths along the whole range poten- 
tially available for A 1500 for a sample of objects, even when 
dust attenuation can be described in the same way for all these 
objects. With these caveats in mind, we determine the mean 
value of ptohe (P) = -1.05 + 0.13 for the lUE sample and 
(P) = -0.50 + 0.11 for the FDF sample (see Tabled}. If the stel- 
lar populations, the dust/stars configuration, and the structure 
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Fig. 9. The slope parameters yi to 74, 712 - yi - 72, and 734 = 73-74 versus f3 for the FDF sample of 2 < z < 2.5 galaxies 
(open circles with eiTor bars) and the lUE sample of local starburst galaxies (crosses). The spectroscopic data are compared to 
combined M05 stellar population synthesis and WGOO radiative transfer models. Details of the models shown are indicated in 
the legend of subfigure (a). The free parameter of the curves is the average radial visual optical depth Ty (see Fig.^. 
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Table 7. UV continuum slope parameters of our sample of 34 FDF galaxies at 2 < z < 2.5 (see Table[0 



No. 


P 


71 


72 


73 


74 


712 


734 



1208 -0.25 ±0.11 -0.72 ±0.12 +0.03 ± 0.08 

1691 -0.29 ±0.11 -0.28 ±0.13 -0.53 ±0.14 

1744 -0.65 ±0.13 +0.68 ±0.14 -0.52 ±0.14 

1991 -1.44 ±0.23 -0.76 ±0.52 -2.24 ± 0.23 

2274 -0.74 ±0.07 -0.17 ±0.05 -0.98 ± 0.05 

2418 -0.21 ±0.09 -0.29 ±0.08 -0.42 ±0.13 

2495 -0.22 ±0.06 -0.03 ± 0.07 -0.56 ± 0.07 

2636 +0.93 ±0.14 +2.37 ±0.15 +0.38 ± 0.09 

3005 -0.14 ±0.09 +0.16 ±0.07 -0.18 ±0.09 

3163 -0.19 ±0.08 +0.44 ±0.07 -0.74 ± 0.07 

3300 -1.36 ±0.08 -1.28 ±0.08 -1.43 ±0.11 

3374 -1.63 ±0.06 -1.48 ±0.07 -1.56 ±0.07 

3688 -0.33 ±0.11 +1.19 ±0.09 -1.88 ±0.10 

3810 -0.33 ±0.06 -0.03 ± 0.04 -0.70 ± 0.09 

3874 +0.01 ± 0.09 +0.94 ± 0.08 - 1 .07 ± 0.09 

3875 -1.52 ±0.11 -1.33 ±0.10 -0.61 ±0.12 
3958 -0.65 ±0.10 -0.18 ±0.11 -1.75 ±0.10 
4795 -0.78 ± 0.08 +0.5 1 ± 0.07 - 1 .04 ± 0.06 
4871 -0.09 ±0.09 +0.27 ± 0.09 -0.61 ±0.09 
4996 +0.03 ±0.07 +1.00 ±0.05 -0.50 ± 0.06 
5058 -1.03 ±0.07 -0.44 ± 0.05 -0.90 ± 0.07 
5135 -0.37 ±0.13 -0.54 ±0.11 -1.51 ±0.15 
5165 -1.06 ±0.09 -0.83 ±0.09 -1.48 ±0.13 
5190 -1.91 ±0.11 -0.79 ±0.07 -1.56 ±0.13 
5227 +0.02 ±0.15 +1.40 ±0.12 -1.32 ±0.13 
6024 -0.99 ±0.06 -0.64 ± 0.06 -0.84 ±0.06 
6372 -0.79 ±0.07 -0.69 ± 0.06 -0.69 ±0.10 
6407 +1.02 ±0.16 +2.69 ±0.17 -0.40 ±0.12 
6934 -0.26 ±0.09 +0.54 ± 0.08 -0.20 ± 0.09 
6947 -0.10 ±0.08 +0.65 ± 0.07 -0.84 ±0.12 
7029 -0.68 ±0.07 -0.35 ±0.06 -0.79 ±0.18 
7078 -0.41 ±0.19 -0.49 ±0.53 -0.45 ±0.15 
7307 -0.56 ±0.11 -0.19 ±0.13 -1.55 ±0.10 
7342 -0.08 ±0.16 -0.42 ±0.15 +0.28 ± 0.20 



-1.21 ±0.09 -1.28 ±0.37 -0.75 ±0.15 +0.06 ±0.39 

-1.02 ±0.16 -0.42 ±0.20 +0.25 ±0.19 -0.60 ± 0.23 

-1.87 ±0.19 -4.30 ±0.38 +1.21 ±0.19 +2.42 ±0.41 

-2.43 ±0.24 -3.72 ±0.79 +1.48 ±0.49 +1.29 ±0.79 

-1.26 ±0.14 -1.68 ±0.43 +0.81 ±0.07 +0.42 ± 0.40 

-1.60 ±0.10 +0.35 ±0.19 +0.13 ±0.16 -1.94 ±0.22 

-2.05 ± 0.07 +0.46 ± 0.69 +0.53 ± 0. 10 -2.5 1 ± 0.69 

-2.83 ±0.15 +1.67 ±0.63 +1.98 ±0.17 -4.49 ±0.71 

-2.34 ±0.10 -0.65 ±0.27 +0.33 ±0.11 -1.69 ±0.29 

-2.34 ±0.09 -0.07 ±0.82 +1.19 ±0.10 -2.27 ± 0.82 

-2.08 ±0.15 -3.27 ±0.66 +0.15 ±0.14 +1.19 ±0.64 

-1.44 ±0.11 -1.23 ±0.27 +0.08 ±0.10 -0.21 ±0.29 

-3.54 ±0.16 -2.96 ±0.59 +3.07 ±0.13 -0.58 ±0.61 

-2.94 ±0.07 +0.30 ±0.24 +0.67 ± 0.09 -3.24 ±0.25 

-4.16 ±0.11 -1.43 ±0.99 +2.01 ±0.12 -2.73 ±1.00 

-2.47 ±0.12 -1.34 ±0.36 -0.72 ±0.16 -1.14 ±0.35 

-2.03 ±0.13 -3.13 ±0.27 +1.56 ±0.15 +1.09 ±0.30 

-1.80 ±0.08 -1.86 ±0.15 +1.55 ±0.09 +0.05 ±0.18 

-2.75 ±0.13 +1.15 ±0.48 +0.89 ±0.12 -3.89 ±0.50 

-1.94 ±0.09 -1.19 ±0.09 +1.50 ±0.08 -0.74 ±0.13 

-1.55 ±0.05 -2.39 ±0.10 +0.45 ± 0.08 +0.84 ±0.11 

-2.75 ±0.18 -2.79 ±0.32 +0.97 ±0.19 +0.04 ±0.34 

-3.03 ±0.17 -2.90 ±0.30 +0.65 ±0.16 -0.14 ±0.37 

-2.69 ±0.24 -1.17 ±0.39 +0.77 ±0.15 -1.53 ±0.38 

-4.85 ±0.18 -0.02 ±0.29 +2.72 ±0.18 -4.83 ± 0.36 

-2.61 ±0.07 -1.33 ±0.35 +0.20 ± 0.08 -1.28 ±0.36 

-1.25 ±0.08 -1.39 ±0.18 -0.01 ±0.13 +0.14 ±0.19 

-0.59 ±0.16 -2.21 ±0.28 +3.09 ±0.21 +1.62 ±0.32 

-3.25 ±0.14 -3.18 ±0.91 +0.75 ±0.12 -0.09 ± 0.92 

-3.45 ±0.12 +0.00 ±0.15 +1.49 ±0.13 -3.46 ±0.19 

-2.99 ±0.14 -0.72 ±0.23 +0.43 ±0.18 -2.26 ± 0.27 

- 1 .8 1 ± 0.42 -4.58 ± 0.77 -0.04 ± 0.54 +2.70 ± 0.98 

-0.34 ±0.17 -2.48 ±0.85 +1.35 ±0.17 +2.19 ±0.86 

-4.05 ±0.20 +0.09 ±0.78 -0.70 ± 0.28 -4.14 ±0.79 



of the dusty ISM are the same for the galaxies of both samples 
(see Sect. 13. 1> . then Fig.|3]implies that the mean attenuation at 
UV wavelengths is larger for the FDF galaxies (see also Noll 
et al. l2004t . In this case, we probably compare z ~ 2 ultralumi- 
nous galaxies (i.e. with Luv <: lO'^ L©) to fainter (by about two 
orders of magnitude) local starburst galaxies (see Heckman et 
al. |T998l Leitherer et al. I2002> . given the gap in observed UV 
luminosities (see Fig.QJ. 

Fig- 13 shows the distribution of the FDF and lUE sample 
galaxies and the loci of the models in several diagrams where 
the individual y parameters are plotted versus fi. Open circles 
and crosses reproduce, respectively, the z ~ 2 UV-luminous 
FDF galaxies and the local lUE starbursts. For both samples, 
the measured parameters have similar uncertainties, though for 
the FDF sample the most uncertain one is 74 due to the bright 
night sky in the corresponding wavelength range. For the lUE 
sample the uncertainties of 73 are particularly high due to the 
low quality of the combined spectra in the overlapping region 
of the lUE short and long wavelength channels (see Sect. 12. 2> . 
Therefore, only error bars for the FDF galaxies are drawn to 
ease the reading of Fig. |51 Dilferent models are reproduced 



synthetically as illustrated in Fig.|5Ja). We discuss the results 
obtained from the analysis of Fig.|5]in the next two sections. 

4.1. The far-UV slope of the extinction law 

Panels (a), (b), and (e) of Fig. |9l show, respectively, 71, 72, and 
712 versus /3 both for the observed spectra and the synthetic 
ones. As the models show, the distribution in these three dia- 
grams characterises the far-UV slope of the extinction law. Not 
surprisingly, both the models and the observed galaxies are dis- 
tributed almost along the one-to-one relation in the /3-yi plane. 
This is due to the large overlap between the wavelength win- 
dows where the two slopes are defined (see Fig. |3- Models 
with different extinction curves split up much better in the ^8-72 
and /3 - 712 planes. In fact, 72 and /3 are defined in wavelength 
regions with only a partial overlap. 

Interestingly, most of the FDF and lUE galaxies are dis- 
tributed in the region delimited by the models with MW or 
SMC dust in each diagram. However, the baricenter of the 
galaxy distribution is close to the locus of the models with 
SMC extinction law, despite the scatter of the distribution is 
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Table 8. Average UV continuum slope parameters of the local lUE sample and different 2 < z < 2.5 PDF samples. The slope 
parameters were averaged unweighted. The uncertainties given are mean errors and were derived from the variance of the sample 
without considering the individual errors. 



Par. 


lUE 


PDF 


PDF 


PDF 


PDF 


PDF 




all 


all 




/3 < -0.4 


-0.4 </3 < 0.1 


-0.4</?<0.1 


-0.4</?<0.1 














734 > -2 


734 < -2 


A' 


24 


34 




16 


16 


8 


8 




-1.05+0.13 


-0.50 + 


0.11 


-1.01+0.11 


-0.18 + 0.03 


-0.23 + 0.04 


-0.12 + 0.04 


71 


-0.51+0.16 


+0.03 + 


0.16 


-0.53 + 0.15 


+0.27 + 0.16 


+0.13 + 0.25 


+0.40 + 0.21 


72 


-0.98 + 0.14 


-0.86 + 


0.11 


-1.15 + 0.13 


-0.67 + 0.14 


-0.65 + 0.24 


-0.69 + 0.17 


73 


-1.63 + 0.26 


-2.33 + 


0.17 


-1.98 + 0.18 


-2.76 + 0.27 


-2.21+0.33 


-3.32 + 0.34 


74 


-1.90 + 0.22 


-1.46 + 


0.26 


-2.34 + 0.30 


-0.73 + 0.33 


-1.51+0.46 


+0.06 + 0.26 


7l2 


+0.46 + 0.21 


+0.88 + 


0.17 


+0.62 + 0.17 


+0.94 + 0.27 


+0.78 + 0.40 


+ 1.10 + 0.36 


734 


+0.27 + 0.32 


-0.87 + 


0.35 


+0.36 + 0.36 


-2.04 + 0.40 


-0.69 + 0.27 


-3.38 + 0.31 




Melallicity: 

° o 0.5 Zq 
Zo 
Age: 

— 50 Myr 

500 Myr 
Extinction: 
□ ■ MW 
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Fig. 8. The combined slope parameter 734 versus /3 for the lUE 
sample of local starburst galaxies. The statistical uncertainties 
of the slope measurements are indicated by eiTor bars. For 
734 these are probably smaller than the systematic uncertain- 
ties (see Section 12. 2> . The spectroscopic data are compared 
to combined |M05 stellar population synthesis and lWGCJUI ra- 
diative transfer models. All M05 models shown (see legend) 
were calculated for a standard Salpeter IMF and continous star- 
formation. All WGOO models given (see legend) were calcu- 
lated for a clumpy SHELL dust/stars configuration. The free 
parameter of the curves is the average radial visual optical 
depth TV (see Fig.|4}. 



large. We conclude that the mixture of dust grains present in 
the UV-luminous FDF galaxies at z ~ 2 produces an extinction 
curve with a far-UV slope that is intermediate between those of 
the MW and SMC extinction curves, if not close to the far-UV 
slope of the SMC extinction curve. 



This result is robust against differences in the stellar popu- 
lations induced by metallicity, age, and star-formation history 
of the real galaxies, for the reasons discussed in Sect. 13.2.21 
Small differences do exist between the distributions of mod- 
els and galaxies reproduced in panels (a), (b), and (e) of Fig.|9l 
however. For instance, yi is sensitive to the maximum age of 
the stellar populations (i.e. the age of the model), while 72 
is sensitive to the metallicity of the stellar populations (see 
Table|5}- Therefore, individual galaxies will shift slightly with 
respect to the loci of models with different dust types in the 
P - ji, P - yi, and fi - yn planes according to the true charac- 
teristic age and metallicity of their stellar populations. 

Furthermore, we note that 26 of 34 UV-luminous FDF 
galaxies at z ~ 2 exhibit > -0.8, which can be achieved 
by the two-phase, clumpy SHELL models with SMC dust for 
Tv > 3 - 5, depending on the properties of the stellar popu- 
lation. In this case the dusty clumps are optically thick at UV 
wavelengths, and the UV photons escaping from the system do 
so after scattering through the diffuse interclump medium that 
is still optically thin (see Sect. 13.1.'^ . The values of the optical 
depths of a typical clump and of the diffuse interclump medium 
depend on the filling factor and density contrast assumed, once 
Tv is fixed (see Witt & Gordon 1996 WGOO I. Nevertheless, 
multiple scattering does play a role in the escape of the UV 
photons from the reddest, dustiest galaxies of the FDF sample 
under study. 

4.2. Evidence for the presence of the UV bump 

Panels (c), (d), and (f) of Fig. |9l show, respectively, 73, 74, and 
734 versus /3 both for the observed spectra and the synthetic 
ones. The parameters 73 and 74 are defined in wavelength win- 
dows corresponding to, respectively, the blue and red wings of 
the 2175 A feature present in the average MW and LMC ex- 
tinction curves. As the models show, the parameters 73 and 734 
are the most sensitive to the presence of the carriers of the UV 
bump as determined from an observed spectrum. Furthermore, 
we recall that 73 and 74 are affected by relatively large eiTors 
for, respectively, the lUE and FDF samples. 

Now, Fig. |9jc) and (f) clearly show that the z ~ 2 UV- 
luminous FDF galaxies populate a region delimited by mod- 
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Fig. 10. The UV-bump indicator 734 as a function of the red- 
dening measure p for our local galaxy sample (cross) and dif- 
ferent subsamples of 2 < z < 2.5 FDF galaxies. The complete 
FDF sample of 34 galaxies (open circle) was divided into a 
blue (y6 < -0.4; 16 objects; filled circle) and a red subsample 
(-0.4 < y6 < 0.1; 16 objects; filled square) excluding two very 
red objects with /3 ~ 1. Moreover, the red sample was further 
fragmented into two subsamples of eight objects each differing 
in 734. The filled lozenge and the filled triangle mark galaxies 
with 734 > -2 and 734 < -2, respectively. The 1 cr error bars 
indicate mean errors. The diagram also shows combined stellar 
population synthesis and radiative transfer models as a function 
of the optical depth of the dust. Details of the models shown are 
indicated in the legend of Fig.|9ja). 

els with either a MW- or an SMC-type extinction curve. 
Furthermore, Fig. |9jf) suggests that 734 becomes more nega- 
tive when j6 increases. This trend means that the 2175 A fea- 
ture characterises the redder, i.e. more attenuated and dustier, 
FDF galaxies. Conversely, the lUE nearby starbursts are almost 
equally distributed across the locus of models with SMC-type 
dust. 

These findings are consistent with those illustrated in the 
previous section, based on 7 parameters defined in distinct 
windows corresponding to the short-wavelength region of the 
UV spectral domain. Hence, we conclude that the redder, UV- 
luminous FDF galaxies at 2 < z < 2.5 contain a mixture of 
dust grains producing an extinction curve that has both a far- 
UV slope and a strength of the UV bump that are intermediate 
between those of the SMC- and MW-type extinction curves. 

4.3. Different dust properties for different galaxy 
subsamples 

Almost all the nearby lUE starbursts exhibit < -0.4, the 
mean value of this parameter being (J3) - -1.05 + 0.13 (see 
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Fig. 11. Comparison of composite spectra of FDF galaxies in 
the range 2.0 < z < 2.5 differing in /3 and 734. Upper spec- 
trum in the upper panel: composite of 16 FDF galaxies with 
/3 < -0.4. Lower spectrum in the upper panel (offset by -0.5 for 
clarity): composite of eight FDF galaxies with -0.4 < yS < 0.1 
and 734 < -2. The lower panel gives the ratio of both com- 
posites. The positions of prominent UV lines and the UV dust 
feature are indicated. 

Table For this reference sample, (734) - +0.27 + 0.32, i.e. 
the evidence for the presence of the UV bump is weak. 

Hence, we divide the z ~ 2 UV-luminous FDF galaxies in 
two subsamples, with either < -0.4 or -0.4 < y6 < 0.1. We 
find that the FDF subsample of 16 objects with /? < -0.4 ex- 
hibits (J3) = -1.01 + 0.11 and (734) = +0.36 + 0.36, in good 
agreement with the reference sample at z ~ 0. Conversely, 
the FDF subsample of 16 objects with -0.4 < /? < 0.1 ex- 
hibits (J3) = -0.18 ± 0.03 and (734) = -2.04 ± 0.40. This con- 
firms that the carriers of the 2175 A dust absorption feature are 
present in red (i.e. with /3 > -0.4) UV-luminous FDF galaxies 
at z ~ 2. Nevertheless, objects with MW-type dust seem not 
to be present among the FDF sample galaxies as well as in the 
lUE sample. 

The presence of the UV-bump can be appreciated directly 
from the composite spectrum of the eight FDF galaxies with 
/3 > -0.4 and 734 < -2 «734) = -3.4), reproduced in Fig.FTTI 
As a reference, the composite spectrum of the 16 FDF galaxies 
with y6 < -0.4 ((734) = 0.4) is also shown. The two spectra are 
normalised in the wavelength range from 1250 to 1500 A. The 
ratio between the "blue" and "red" composite spectra exhibits a 
broad feature centred at 2166 ± 15 A, consistent with the centre 
of the UV bump. Its slope is negative owing to the different 
values of /3 for the two subsamples. 

Having established the presence of a UV bump in the ex- 
tinction curves of the red UV-luminous FDF galaxies at z ~ 2, 
can we also estimate its strength? Unfortunately, there are no 
I WGOOi models for extinction curves different from those typi- 
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Fig. 12. The UV-bump indicator 734 as a function of the red- 
dening measure yS for our local (crosses) and 2 < z < 2.5 
FDF galaxies (circles). The diagram also shows combined stel- 
lar population synthesis and dust attenuation models as a func- 
tion of the optical depth of the dust. The dust attenuation mod- 
els plotted are the Calzetti et al. (2000 1 law (long dashed lines 
and ckcles) and the .WGOO. models for the clumpy SHELL con- 
figuration with SMC (solid lines) and MW dust (dash-dotted 
lines). Furthermore, the extinction laws of the SMC bar (dot- 
ted lines and lozenges), the LMC 2 supershell (dashed lines and 
lozenges), the LMC average (dashed lines and squares), and the 
MW (dotted lines and squares) are indicated fCCM'; Gordon et 
al. .2003j . The symbols are plotted in intervals of AEb-v = 0.1. 
In all cases a lM05l model with Salpeter IMF, continuous star- 
formation, an age of 200 Myr, and solar metallicity has been 
taken. 



cal either of the diffuse ISM of the Milky Way or of the SMC 
bar. Therefore, for illustrational purposes only, we show the 
distribution in the j3 - 734 plane of models where dust attenua- 
tion is described by a non-scattering homogeneous dust screen 
for different empirical extinction laws as a function of Eb-v 
(see Calzetti et al. 120001 Gordon et al. 2()0Jl. In particular, 
we consider the extinction functions determined for the SMC 
bar, the LMC 2 supershell region, the average LMC (Gordon et 
al. '20031 and the average MW fCCM), that are displayed in 
Fig. As a reference, we take the attenuation functions pro- 
duced by the two-phase, clumpy SHELL models with MW- 
or SMC-type dust ofWGOO' and the Calzetti et al. ('2000^ at- 
tenuation law. This comparison holds only to the purpose of 
comparing slopes of synthetic spectra (see Sect. l3.1.2t . Fig. 1121 
reproduces the loci of the different models and the distribution 
of the FDF and lUE sample galaxies. It is evident that a mix- 
ture of dust grains with properties ranging from those typical 
of the SMC bar to those typical of the LMC 2 supershell region 
can account for the distribution of the data in the j0 - 734 plane. 



Fig. 13. The EW of the metallicity indicator C IV as a func- 
tion of the UV-bump strength indicator 734 for our local galaxy 
sample and different subsamples of 2 < z < 2.5 FDF galaxies 
(see legend and caption of Fig.llO>. The EWs of C IV were de- 
rived from composite spectra of the different samples. The 1 cr 
eiTor bars include the variance between the individual spectra, 
the errors of the continuum level definition, and the statistical 
uncertainties of the line strengths. 

We recall that the LMC 2 supershell region is located near the 
extremely vigorous 30 Dor star-formation region (Misselt et al. 
1999 1. Conversely, extinction curves typical of more quiescent 
star-formation activities like those of the average LMC and 
MW seem not to describe the distribution of the data in the 
previous plane. This is consistent with the nature of actively 
star-forming galaxies of the FDF and lUE sample galaxies (see 
Sect.O- 

Comparing the 712 values for the different subsamples de- 
fined above (see Table |8}, we do not find a significant relation 
between the far-UV slopes of the extinction function as indi- 
cated by 712 and the bump strength. The moderate increase of 
712 with /3 is consistent with the predictions from the IWGOOl 
models for SMC-type dust, for which the mean values are in 
better agreement than for the MW-type models (see Sect. l4.lt . 

4.4. Presence of the UV bump and strengths of stellar 
and interstellar lines 

Hereafter we investigate whether the presence of the UV bump 
is linked either to the (stellar) metallicity or to properties of 
the ISM of the system, as traced by the equivalent widths of 
the absorption lines listed in Tab. |3] We group the objects in 
subsamples, as in Sect. 14.31 

Fi g . 1 1 3 Ireproduce s 734 versus Wciy. Wciv is a good metal- 
licity indicator in starburst galaxies, the CIV feature being 
mainly produced by hot star winds with a strength that de- 
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Fig. 14. The average strength of six prominent interstellar low- 
ionisation absorption lines (see Noll et al.5004 1 as a function of 
the UV-bump strength indicator 734 for our local galaxy sample 
and different subsamples of 2 < z < 2.5 FDF galaxies (see 
legend and caption of Fig.llO>. The EWs of the low-ionisation 
lines were derived in the same way as reported in Fig. ll3lfor 
CIV. 



pends on the chemical composition (Walborn et al. 1995 
Heckman et al. [T998l Leitherer et al. l20OTl Mehlert et al. 2002 1. 
Interestingly, Wciv is smaller for the subsample of blue UV- 
luminous FDF galaxies at z ~ 2, as well as for the lUE sample 
of nearby starbursts, than for the red subsample of FDF galax- 
ies. However, red FDF galaxies selected according to their val- 
ues of 734 do not exhibit significant differences in Wcw- Hence, 
the metallicity of the OB stars (and, thus, of the ISM where 
these stars originated) does not primarily determine the pres- 
ence of the carrier of the 2175 A dust absorption feature among 
the z ~ 2 UV-luminous FDF galaxies. This is consistent with 
studies of extinction curves in the local universe (see Sect.[0. 
On the other hand, metallicity is higher (by more than 2cr) 
for FDF objects affected by a larger reddening at rest-frame 
UV wavelengths (i.e. with a larger /?). This result holds for the 
whole FDF spectroscopic sample (Noll et al. l2004t . 

Finally, the parameter 734 seems to depend also on the av- 
erage EW of the six most prominent low-ionisation absorp- 
tion lines that fall within the rest-frame UV domain mapped 
by FORS (VKlis, see Fig. ll4> . In particular, values of 734 close 
to zero are associated with low values of Wus- They charac- 
terise the lUE local starbursts and the blue z ~ 2 UV-luminous 
FDF galaxies. On average, the red z ~ 2 UV-luminous FDF 
galaxies exhibit larger values of Wus (by more than 1.5 cr). 
This difference is essentially caused by the red FDF objects 
with 734 < -2, having larger values of Wus at almost the 3 cr 
level. 




-1 -0.5 
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Fig. 15. The logarithm of the total stellar mass (Drory et al. 
12005^ in solar units (upper panel) and the half-light radius 
(Pannella et al. 2005 1 in kpc (lower panel) as a function of the 
reddening measure /3 for different subsamples of 2 < z < 2.5 
FDF galaxies (see legend and caption of Fig.llO>. The 1 cr er- 
ror bars are based on the variance in mass and radius between 
different galaxies. 

This behaviour can be interpreted as follows. The low- 
ionisation absorption lines trace the spatial and velocity dis- 
tribution of the neutral interstellar gas illuminated by the rest- 
frame UV continuum along the observer's line of sight. In par- 
ticular, Wlis is a proxy for the neutral clouds' covering fraction 
(Shapley et al. 20Q3|. Larger values of Wus imply a larger ob- 
scuration at rest-frame UV wavelengths, if the spatial distribu- 
tions of neutral gas and dust are correlated. Hence, they imply 
more negative values of 734, if the extinction curve exhibits a 
UV bump, as shown by Fig. ^| for the dusty media modelled 
bv WGOO. 

5. Discussion 

Dust grains play a fundamental role in the formation and evo- 
lution of galaxies. In fact, they are both a byproduct of stel- 
lar evolution (in particular, a sink of metals in the gas phase), 
as well as a site for the efficient formation of H2 molecules 
(e.g. Cazaux & Tielens '2002' '2004 1. Furthermore, dust grains 
absorb stellar light and re-emit it at mid-IR-mm wavelengths, 
controlling the energy balance in the ISM and protostellar gas 
clouds on one hand, and complicating the interpretation of ob- 
served SEDs of any dusty system on the other. Finally, they can 
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offer the surface that is needed for momentum-driven winds on 
both stellar and galactic scales. Those effects of dust already 
seem to take place when the metallicity of the environment is 
very low (i.e. ~ 0.01 Zq, Hirashita & Ferrara 2002 Morgan 
& Edmunds I2003> . For all those reasons, the first dust enrich- 
ment in the course of galaxy evolution may drive an enhanced, 
high star-formation rate, shaping the cosmic SFR history at the 
highest redshifts. 

The first source of dust in the universe is Type 11 supernovae 
or pair instability supernovae, beginning to occur a few Myr af- 
ter the onset of star-formation. After a few 10^ Myr evolved in- 
termediate/low mass stars start to contribute to the dust enrich- 
ment. Hence, for the individual galaxy, the relative importance 
of the two sources of dust depends on time and SFR history 
(see, e.g., Dwek 19981. 

The extinction curves expected in young (i.e. less than 
1 Gyr old) galaxies are quite sensitive to internal metal mix- 
ing in supernovae (Hirashita et al. l2005> . In general, these ex- 
tinction curves do not exhibit a 2175 A absorption feature (as 
the standard MW extinction curve instead) because the con- 
tribution of carbon grains is small (for dust yields, see Todini 
& Ferrara .2001. Nozawa et al. 2003, Schneider et al. .2004 ). 
Furthermore, they are different from the MW, LMC, and SMC 
extinction curves, especially at rest-frame UV wavelengths 
(Hirashita et al. l2005> . This is no surprise, since the ages of the 
MW, LMC, and SMC are much longer than 1 Gyr and, thus, 
at least a change in the dust production mechanism has taken 
place. Furthermore, the dusty gas outflows of intermediate/low 
mass stars are carbon rich. 

At present, there are only few observational studies on the 
extinction curves in systems at very high redshift. The extinc- 
tion properties of broad absorption line quasars at z > 4.9 
seem to be different from those of analogous objects at z < 4 
(Maiolino et al. .2004a, i. In particular, the extinction curve of 
SDSS J 1048-1-4637 at z = 6.2 has been shown to be fitted by 
dust models of Type II supernovae rather than those of pair in- 
stability supernovae (Maiolino et al. .2004b., i. For ages of the 
universe much longer than 1 Gyr (i.e. for z « 5), the known 
extinction curves provide controversial conclusions as for the 
presence of the 2175 A absorption feature (see Sect. This 
may be a consequence of at least a difference in the "character- 
istic" SFR history of different classes of objects. 

In particular, the absence of a pronounced UV bump seems 
to characterise the average extinction curve in the large, well 
studied photometric sample of LBGs at z ~ 3 observed by 
Steidel et al. (2003) (Vijh et al. .2003 ). Interestingly, the UV 
bump is weak if not absent also in starburst galaxies of the lo- 
cal universe, whatever their metallicity (Calzetti et al. '1994, 
Gordon et al. 1997) . Is this so for all UV-luminous galaxies 
down to z = then? 

In order to answer to this question, we have investigated the 
extinction properties in 34 UV-luminous galaxies at 2 < z < 2.5 
from the FDF spectroscopic survey (see Sect. I2.lt . Using a 
new parametric description of the rest-frame UV spectra (see 
Sect. 13. 2> we have found that most of the z ~ 2 UV-luminous 
galaxies are characterised by an extinction curve that is inter- 
mediate between those of the SMC and the MW (see Sect.l^. 
This means that the characteristic extinction curve of these ob- 



jects exhibits a weaker 2175 A feature but a steeper far-UV rise 
with respect to the average extinction curve for the diffuse ISM 
of the MW. 

The bulk of the mixture of dust grains (regarding size 
and chemical composition) present in all the 34 UV-luminous 
galaxies at 2 < z < 2.5 is different from the mixture of dust 
grains producing the standard MW extinction curve. This is 
true especially for those galaxies that have probably a low dust 
content (and, thus, a low value of as for our models). On the 
other hand, SMC-like dust may characterise up to 30% of these 
galaxies. In conclusion, the bulk properties of dust in z ~ 2 UV- 
luminous galaxies seem to range between those in the SMC and 
LMC. 

The size of the z ~ 2 UV-luminous galaxy sample allows 
sub-sampling. In general, the 18 "red" (i.e. with fi > -0.4) UV- 
luminous galaxies at 2 < z < 2.5 exhibit a stronger evidence of 
the presence of the UV bump. The presence of the UV-bump in 
the eight UV-luminous galaxies at z ~ 2 with -0.4 < < 0.1 
and < -2 ((734) = -3.4) can be appreciated directly from 
the ratio of the "blue" and "red" FDF composite spectra repro- 
duced in Fig. ^2 This result is at variance with that obtained 
for the nearby starburst sample, that appears to be characterised 
by SMC-like dust whatever fi, consistent with Gordon et al. 
( 1997), though it mostly exhibits (3 < -0.4. Our interpretation 
of this finding is that the caiTier of the 2175 A feature is more 
abundant in redder and, thus, dustier objects in the z ~ 2 UV- 
luminous galaxy sample. 

On the other hand, z ~ 2 UV-luminous galaxies with (3 ~ 
exhibit a large scatter in both the parameters jn and 734 (see 
Fig. 1121 and Table Now, 712 is particularly sensitive to the 
far-UV slope of the extinction curve, while 734 is particularly 
sensitive to the presence (and strength) of the UV bump. We 
conclude that there is no evidence of a correlation between the 
strength of the far-UV slope and the strength of the UV bump, 
for the extinction curves of the FDF galaxies with evidence of 
the UV bump. This may suggest that different environments 
(in terms of properties of the ISM and of the general radiation 
field) and/or different degrees of dust processing characterise 
different z ~ 2 UV-luminous galaxies, if the lessons derived 
from the local universe apply to the universe at z ~ 2 (see 
Sect.Hl). 

The available data allow us to check the first suggestion. 
We find that the FDF sample galaxies with /? ~ are charac- 
terised by relatively high metallicities of Z ~ 0.5 - 1 Zq (see 
Mehlert et al. 120021 l2005l . However, flie metalUcity indicator 
Wciv does not show any significant correlation with 734 (see 
Fig. I13> . Hence, a relatively high metallicity (i.e. close to the 
solar value) seems to be a prerequisite for the occuiTence of the 
UV bump. However, it does not determine the strength of the 
2175 A feature. 

Furthermore, the z ~ 2 UV-luminous galaxies with ev- 
idence of a relatively strong UV bump are among the most 
massive and largest galaxies at this redshift. This is shown in 
Fig- El where we plot the estimated total stellar mass^ (upper 

^ The total stellar mass was estimated by Drory et al. <2005> fol- 
lowing Brinchmann & Ellis ( 2000 ). In particular, Drory et al. <2005> 
derived stellar masses after applying a large set of two-component 
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panel) and the rest-frame UV effective radius^ (lower panel) 
versus yS. Nevertheless, there is no coiTelation between stel- 
lar mass or size of a galaxy and the strength of the UV bump. 
In conclusion, an extinction curve exhibiting the UV bump is 
likely to be associated with the most massive, largest, and, pos- 
sibly, more evolved systems among the z ~ 2 UV-luminous 
ones. 

Finally, we explore the possibility that the topology of the 
dusty ISM produces the large scatter of 734 observed for the 
red UV-luminous galaxies at 2 < z < 2.5. This hypothesis is 
made appealing by the fact that the red UV-luminous galaxies 
with < -2 exhibit a significantly larger Wus than the aver- 
age z ~ 2 UV-luminous galaxy and the average local starburst 
galaxy. We recall that Wlis traces the covering fraction of the 
young stellar population by (dusty) neutral gas clouds in the di- 
rection to the outside observer A larger covering fraction may 
be due to inclination effects, for a disk geometry (e.g. Pierini 
et al. l2004t . For the z ~ 2 FDF objects with -0.4 < yS < 0.1 
we find small but insignificant differences in the average el- 
lipticities e ^ for the 734 > -2 and 734 < -2 subsamples (see 
TableEJ, indicating (e) = 0.41 + 0.07 and 0.55 + 0.08, respec- 
tively. Hence, larger values of Wlis possibly imply a better self- 
shielding of the dust grains that populate the intervening gas 
clouds, despite inclination effects may be not negligible. The 
carriers of the 2175 A feature would then be protected from 
strong and hard radiation fields by a screen of other (more ro- 
bust) dust grains (e.g. Gordon et al. 2003- In this case, light 
emitted by young stars and scattered by the embedding dusty 
gas clouds plus light emitted by stars living in relatively qui- 
escent environments (e.g. the diffuse ISM) would dominate the 
UV spectra. On the other hand, smaller values of Wlis may im- 
ply a larger dust-clearing effect by superwinds (e.g. Strickland 
& Stevens I2000t in the direction to the outside observer. The 
effect of superwinds is probably larger in less evolved systems. 
In this case, the UV spectra would be dominated by light pro- 
duced in star-formation regions, characterised by strong radia- 
tion fields and shocks. In these harsh environments, probably 
the caiTiers of the 2175 A feature would have a larger probabil- 
ity to be destroyed, if ever produced. However, the local uni- 
verse shows that the latter conclusion is not so straightforward 
(see Sect-HJ- 



models to the whole FDF photometric data (see Heidt et al. 12003 > . 
Each model consisted of a combination of two composite stellar pop- 
ulations, as given by a young (lOOMyr) burst with a constant SFR 
plus an older burst (> 500 Myr) with an exponentially declining SFR. 
Both stellar components were attenuated independently according to 
the "Calzetti law" (Calzetti et al. 2000). 

* Effective radii and ellipticities were estimated by Pannella et al. 
<2005il . who determined the morphology of the whole FDF galaxies 
by applying the publicly available two-dimensional surface brightness 
profile fitting packages G1M2D (Simard et al. [19991 and GALFIT 
(Peng et al. 2002 1 to publicly available Hubble ACS images of the 
FDF in the F814W broad-band filter. The effective wavelength of this 
filter corresponds to rest-frame wavelengths of about 2500 A at z ~ 2, 
as for the sample under investigation. 

^ e = -\/l - b-/a^, a = major axis, b = minor axis; for the origin of 
the data see Pannella et al. 12005 j and footnote 6 



6. Conclusions 

We have investigated the properties of the extinction curve in 
the rest-frame UV for a sample of 34 UV-luminous galaxies 
at redshift 2 < z < 2.5 selected from the FORS Deep Field 
spectroscopic survey (Noll et al. 2004 j . Our analysis makes 
use of a new parametric description of the rest-frame UV spec- 
tra of the observed objects and of synthetic spectra. The latter 
spectra are computed with models combining stellar popula- 
tion evolutionary synthesis (Maraston 2005, M05 1 and radia- 
tive transfer of the stellar and scattered radiation through dif- 
ferent dusty interstellar media for the SHELL dust/stars config- 
uration, in a spherically symmetric geometry (Witt & Gordon 
2000, WGOOt. 

Firstly, we have established how much the new parameters 
(seven in total) depend on the properties of the stellar popula- 
tions or of dust attenuation assumed for galaxy models aimed 
at reproducing the rest-frame UV spectra of star-forming galax- 
ies. 

Secondly, the robustness of our method has been 
tested against a reference sample of 24 nearby starbursts. 
Comfortably, we reproduce the result that local starburst galax- 
ies possibly contain dust with an extinction curve lacking the 
2175 A absorption feature (or "UV bump"), like the Small 
Magellanic Cloud (SMC) but at variance with the Milky Way 
(MW), and a steep far-UV rise, intermediate between those in 
the SMC and the MW extinction curves (Calzetti et al. 119941 
Gordon etal. [T997l . 

The application of the same method to the sample of z ~ 2 
UV-luminous galaxies leads to the following results. 

- The z ~ 2 UV-luminous galaxies at 2 < z < 2.5 with a 
relatively low amount of reddening at rest-frame UV wave- 
lengths (i.e. with a far-UV slope /? < -0.4) contain a mix- 
ture of dust grains that is similar to that held to be present in 
local starburst galaxies. This analogy holds despite the fact 
that the former galaxies are about one order of magnitude 
more luminous than the latter 

- For the highly reddened (i.e. with -0.4 <y6<0.1)z~2 
UV-luminous galaxies, the extinction curves range between 
those typical of the SMC and Large Magellanic Cloud 
(LMC). Therefore, part of these objects exhibits UV bumps 
as pronounced as in the LMC. 

- The amount of reddening at rest-frame UV wavelengths is 
larger in more metal-rich objects, the equivalent width of 
the CIV doublet (Wciv) being used as an indicator of the 
total metallicity. On the other hand, objects with large val- 
ues of Wciv (i-C- metal rich) can exhibit both SMC- and 
LMC-type dust, despite the 2175 A absorption feature is 
associated with the presence of carbonaceous dust grains. 

- The presence of the UV bump seems to be associated with 
a large average equivalent width of the six most prominent, 
rest-frame UV, low-ionisation absorption lines falling in the 
FORS spectra. The average equivalent width of these satu- 
rated lines offers a proxy for the covering fraction of young 
massive stars by neutral gas clouds (containing dust grains), 
and, thus, for the ISM topology. 

- Furthermore, we find that the most opaque and, thus (for 
our models), dustiest UV-luminous FDF galaxies at z ~ 2 
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tend to be among the most metal rich (i.e. with an average 
metallicity (Z) ~ 0.5 - 1 Z©), most massive (i.e. with an 
average stellar mass (Mstars) ~ 6 x 10''^ Mq), and largest 
(i.e. with an average rest-frame UV effective radius (/?eif ) ~ 
4kpc) systems. 

We interpret these results as the evidence for a difference in 
the properties of the dusty ISM among the most evolved UV- 
luminous, massive galaxies at z ~ 2. As for the local universe, 
the diversity in the bulk properties of the dust grains from one 
galaxy to another is possibly caused by the sum of the local 
effects due to stellar radiation field, shocks, ISM topology, and 
history of chemical enrichment. In particular, the degree of self- 
shielding by dust may play an important role for the production 
and/or survival of the small carbonaceous grains that are held 
to be the carriers of the UV bump already at z ~ 2. 
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